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ABSTRACT 


This thesis describes a procedure for computing a lower confidence 
limit on the reliability of a quasi-coherent complex system using test data 
On its components. The failure times of the components are assumed to 
have either exponential or Weibull distributions with unknown parame- 
ters. The accuracy of this procedure is evaluated using computer simu- 
lation for various system structures and sets of parameter values for the 
assumed distributions. This thesis is an extension of a thesis by Kah Chee 
Yee in that it uses a different equation for the estimate of the shape pa- 
rameter in the Weibull distribution than Yee used, and it evaluates the 


procedure for a larger collection of system structures. 
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I. INTRODUCTION 


This thesis modifies an existing approximate lower confidence limit 
procedure for the reliability of complex systems developed by Yee 
| Ref. 7]. The purpose of the modification is to make the procedure more 
accurate for more configurations of quasi-coherent complex systems and 
easier to use computationally. A system is defined to be quasi-coherent 
if an increase in reliability of any one of its components does not cause 
a decrease in the system reliability. The components of a quasi-coherent 
system do not need to be Statistically independent. However, throughout 
this thesis it is assumed that all components are statistically independent 
and the probability distributions of their failure times are either expo- 
nential or Weibull . 

The approximate system reliability lower confidence limit procedure 
developed by Yee is quite accurate for series systems. In this thesis the 
procedure for estimating the shape parameter, f, differs from the maxi- 
mum likelihood method used by Yee. The procedure developed by 
Varadan [ Ref. 6] is used to estimate the shape parameter, f, in the 
Weibull distribution in this thesis. 

In addition to modifying the lower confidence interval estimation 
procedure developed by Yee, more complex structures are examined here 
than in Yee’s thesis. The computer programs developed by Yee were 
modified to examine these new structures. Also, a computer program was 
developed that can be used to compute the lower confidence limit for the 


reliability of a complex system using these procedures. The system is de- 


scribed by the user in response to queries by the program. Also, the test 


data set is entered by the user in response to queries. 


II. ESTIMATES FOR THE SHAPE PARAMETER OF A WEIBULL 
DISTRIBUTION 


If the time to failure, X, of a device has a Weibull distribution with 
scale and shape parameters @ and f respectively, then its probability 


density function 1s given by 
a(x; 8, B) = (1/6) Bx? exp{—(x/)"}, x>0 (2.1) 


This property is stated more briefly by the phrase X is WEI( 6, f ). The 


cumulative distribution function, CDF for X 1s 

G(x; 0, B)=1 — exp{ —(x/6)°}, x>0 (2.2) 
and the survival function 1s 

G(x; 6, B)=1 — G(x; 9, B) (2x3) 


It is well known, see [ Ref. 6], that the random variable Y=In YX has the 


extreme value distribution with density function f and CDF F given by 


fyi ps0) = eM? exp{ cP}, 00 <p < 00 
aaa Tie-aiore (2.4) 
ag>Q 
Fy; vw, o) = 1 — exp{ —eV-HIo, (255) 


where o = 1/8, 4=1n @ and In is the natural logarithm. In this case we 
write Y is EV( oa, p ). 


Suppose X),.4),...,4, are independent random variables with a common 
WE]I( @, B ) distribution and let X,)< Xy) < ... << Xj, be the first r ordered 
set of the original set of n variables. Then the following equations are 
solved to obtain the maximum likelihood estimates (MLE) f and B for @ 
and B, see [ Ref. 1] : 


- yf p 
Fo X(j) In X(;) a (11 = r)X(,) In X(r) ly 
—<—> > =, * eee (2.6) 


3 p72 
Fx; X(j) + (n —") Xr) / 


(p= (2.7) 
2 X(j) + (n ae 


Closed form expressions for do not exist, and iterative procedures are 
used to solve for B and 8. Computer programs are readily available to 
compute B and 0. They require an original estimate f, to start an iterative 
process, and if it is not chosen carefully the iteration will not yield an 
accurate estimate for f. The MLE estimate, B, is biased. Bain and Lee 
[ Ref. 5] have an excellent treatment of the Weibull distribution and 
properties of the MLE estimates 9 and B. Their discussion includes a table 
of factors (page 200) which can be used to construct nearly unbiased es- 
timates for fp. 

Balakrishnan and Varadan [ Ref. 6 | have derived a method for esti- 
mating f and @ which does not require computer iterations. Their method 
will be used in this thesis to modify Yee’s existing procedure that derives 
lower confidence limits for the reliability of complex systems that have 
some components whose failure times have a Weibull distribution. A 


summary of their results follows. 


Let Yost Ygaays--+> Vines) be a sequence of the order statistics for a ran- 
dom sample of size n variables each with density function given by (2.4). 
Then 


Yo+l) = Yr+2) Sos Vin=s) (2.8) 


If we let X, = (Y)—yp)/o, then the likelihood function based on the 


censored sample Is : 


a A rcs ces 
— ang Laker alee Xa). AM (2.9) 
where, 
A =e 8 
F(x) = 1—exp(—e’) 


F(x) = 1- F(x) 
f(x) = e*exp(—e*) 


The likelihood equations for uw and oa are: 


PiniL p  AXr41) Mee SXa-9) LEO 5 
eae th) I, Fey Ho 210 


=——|r 
ou : f (Xt)  F( (Xing) “la 
SInL MXG-+1) 
- Se) ae 2.1] 
da g lA + 7X41) FUX(y41)) . 





IX n—s)) 4 . F(X) 
+ 2 to Fay I=? 


=r+l1 


Equations (2.10) and (2.11) do not admit explicit solutions. Let, L* denote 
the approximated likelihood function. Using a Taylor series expansion for 
In L to obtain In L*, the two equations can be approximated by : 


x 
SP SR = Fle — 8X41) - (2.12) 


n—S 
s(l — Oy—s Bn—sX(n—s)) er (a; a BX] tl 


dInL _ din L” 


50 ~ ve =— Slt Xan — OXyay) - (2906) 


um—S 
SX(y—5\(1 O_o 6 T Bn—sX(n—s)) a Xi; — BX] ==. 
See Balakrishnan and Varadan [ Ref. 6 p-147] 


Solving (2.12) and (2.13) for yw and o gives their approximate maxi- 


mum likelihood estimates as follows : 





N= B-8C (2.14) 
¢={—D+(D? + 4AB)4/24 (2.15) 
wirere , 

pple. ie 2.16) 
gq, = 1-p:; (2.17) 
C. Seli+ Ing —inG— ing) (2.18) 








= TG; (2.19) 
6 = Nad) In {] ein } (2.20) 
a Pe | r+] [ten r+] ; 
ee ie in ing..)} + (2.21) 
= Proj G+ Gr+] 
Gr+1 2 

I niggiey elena ewy) 
Pr+] 
A=n-r-s Or?) 
: m—S 

B= {ro Yr+1) Sts Ss) Yn—s) oe Lace lm (2.23) 
n—S 

Ce {ry—sti-a,_)+ 2% a; }/m (2.24) 
i=r+] 

fa 70 + sp. +2 B; (25) 

N—S 
D = ry Yirg1) —B) — 1 = Ay N¥n-s) —B) +E ak Y(y —B) (2.26) 
a 2 a> 2 
E= ro(Y (41) —B)+ SP pes’ (ns) — B) coe BAY — B) Ce) 


The approximate MLE estimates of u and o as defined in equations (2.14) 
and (2.15) are biased estimators. Balakrishnan and Varadan 
[ Ref.6 p - 149 | provide a table of constants that are called SBIAS(n,r,s). 


They are a function of the number n of test items placed on test, the 


number r where the observations began and the number s of successive 
components observed. In this thesis the parameter r will always be zero, 
meaning that observations always start from the smallest failure time. 
Using their table an approximate unbiased estimate o* for o is (See Ap- 
pendix A) 


A 


AX 0 
=> ——— 22 

: SBIAS(n,0,s) (2a 
The inverse of o is B ,which will be a biased estimator for f . Bain 

[ Ref. 5 p-220] provides a table of constants B(n), which depend on the 


number of test items, n, such that Be = B {B(n)} is nearly unbiased for B. 


Ill. DESCRIPTION OF THE LOWER CONFIDENCE LIMIT 
PROCEDURE 


The lower confidence limit procedure developed in this thesis is a 
modification of the procedure in the thesis written by Yee , lipRetaaly: 
The procedure in this thesis uses different estimators for the parameters 
B and @ in the Weibull distribution. The method of estimation presented 
by Mann and others [ Ref. | ] and used by Yee employs the maximum 
likelihood estimates for # and @ which require computer iteration and a 
reasonably good guess for an initial value of f to perform the iteration. 
The estimators developed by Balakrishnan and Varadan [ Ref. 6] is used 
in this thesis. They provide an alternative estimation procedure which 
does not require computer iteration methods to compute the estimates. 

In this thesis we consider systems that are made up of k independent 
component subsystems. Systems undergo missions of some duration, say 
t. During a mission the system components are subject to periods of ac- 
tivity and inactivity. During active periods component i is subject to 
failure, with a hazard rate A(t). During inactive periods the hazard rates 
are all zero. A component successfully completes a mission if its total 
operating time in L0,1] exceeds some design parameter 1#(t). A system 
completes its mission if sufficiently many of its subsystems do; the system 


reliability is, as usual, a function of the structure of the system. 


A. SERIES SYSTEMS _ 
1. Interval Estimation Procedure for Exponential Failure Times 
Suppose a series system has k components whose failure times are 


Statistically independent. Suppose the distribution of the failure time of 


component i is exponential with failure rate J;. Then the system reliability 


R, can be written as a function of /, and ¢, 1 = 1,2,...,k as follows : 
k 
R,{t) = exp( 74; d it; ) ie 1) 


where t;, = t(t) 1s the time component 1 operates when the system operates 
for time t. Using the relationship 7,=4/,/2,,, for i= 1,2,...,k where 4,, is 


the failure rate of any one of the k components, equation (3.1) becomes 
k 
R(t) = exp( — Am™ rt; ) (OY) 


If the values of the 7, are known and if am, | is an upper 100(1 — a) % 
confidence limit for /,,, the corresponding approximate lower confidence 


limit for R(t) would be : 
A A k 
Rt) 1 (a) = exp( a Am,U(a) 2 rt; ) (3%) 


The equation for /,,y,, depends on the plan for testing the components. 
The following case is considered in this thesis : 
If n, items of component type i are tested until f failures occur, 1 


= 1,2,...k and if Xiq), Xiq)s--.X) denote these ordered /; failure times then 


2 
A Xa,2F 
AUG) =n (3.4) 
a r;T; 
[j= 


where 7; denotes the total test time accumulated on all n, items of type 


Si k 
ij ie, T=(4—-fXg + LU Xm F = Lf and y2,- is the 1001 — a)th 
a f=] 


10 


percentile point of a Chi-square distribution with 2F degrees of freedom. 
See Bain and Engelhardt [.Ref. 3 ]. 
Values of the 7; are assumed to be unknown in this thesis. They 


will be estimated by 7,, a nearly unbiased estimator for r,, defined by 


A 
A A; 


= — (3.5) 


i 
“mm 





<> 


where 1, = (f; -1)/T, which is an unbiased estimator for 2, when f/> | 
LRef. 7 p. 39]. If Lf, were unbiased for 1//,, and if aM and oe are inde- 
pendent, then 7, would be an unbiased estimator for r,. Replacing m with 
Lf, —1) in equation (3.5) will yield an unbiased estimator fr, for 1. 
Multiplying by this constant f,,/(,, —1) is nullified by a cancellation with 
the same constant in the final equation for the system reliability lower 
confidence limit, so equation (3.5) is used to estimate 7. Using the esti- 


mator 7, for r,, equation (3.4) becomes 


2 
A Xo.2F 
4m, U(a) = (3.6) 
22 rT; 
p= 


It is important to note that the index m denotes the component for which 
ii = (f;—1)/T, is the largest among all the components in the system. The 
corresponding equation for the 100(1 — «) % lower confidence limit on 


the reliability of the series system becomes 


A A k A 
RA) (a) = exp( = Am,U(a).~, rt; ) Gwe 


2. Interval Estimation Procedure for Weibull Failure Times 
Suppose a series system has k statistically independent components 
and suppose the failure times of all k components have Weibull distrib- 
utions. Let the time to failure, X, , of components of type 1 have density 


function 
Bea bea B; 
Fit, = hi B it; exp{ — (A t;) } ; l; > 0 (3.8) 
Then the system reliability R,(t) can be written as 


k 
RQ) = 1 expl — Git") (3.9) 
k 
= exp| ro jPigP 
— 


is 
= exp[ — Ane rth | 
j= 
where 7% = Af, A* is any one of the 7* , 1=1,2,...,.k and 7 =A*/A~ “Tite 
B; values are known, then Xf: will have an exponential distribution with 
constant failure rate 1% and the procedures described in Section | can be 


used to obtain RD) 0) with 7, defined by 
i (n, — fx?! + ; xb (3.10) 
i { UA i(f;) a i) ; 


Suppose f, is unknown and Xq),Xjq)).-,Ayy) are the ordered failure times 
under failure truncated testing for component 1 in the system. The maxi- 
mum likelihood estimates f, for B, and 4, for 4, are the solutions for B, 


and /,; in equations : 


B; ly 
~~ = SS ln Xq 3.11) 


B; jijal 
Xih + 4 A 
rs ; nme (ori2) 
Xi + (4; PX 
where X;,, = Ay) under failure truncation testing. See Mann and others 


[ Ref. 1 p. 189-191]. 

Equations (3.11) and (3.12), can not be solved in closed form. 
They may be solved using iterative computer methods, but in this thesis 
another method is used which does not need computer iterations to ob- 
tain B, and yy . This method was described in Chapter IJ. It is an ap- 


proximation of the maximum likelihood estimate. Let, 
(Balls) 


In this thesis the distribution of 7,, 


y 1S approximated by the exponential 


distribution with failure rate Ae: = Ax. 
The procedures for obtaining the lower confidence limit for system 


reliability are similar to those in Section A. Define 


~e> 
may x< 


r= (3.14) 





~~ > 


x x 


where, 





eam = 1; ’ i= rey (Sila 
rf 
jad 
and 
2% = max 2%. 
l 
Then equation (3.14) becomes 
Ron Jb 
r= 1i(F>) (3.16) 
Hr 


The approximate 100(1 — «) % upper confidence limit for 4% is : 


2 
By. Xo .2F* 
0) = (3.17) 
25 AT, 
I 
where, 
k 
FX = 2 Si (3.18) 
i 


The corresponding approximate 100(1 — a) % lower confidence limit for 


the reliability R(t) of the series system is given by : 
: an 
RO) 1 (a) = exp| — 2 n,U(a)™, it; a (3.19) 


The accuracies of these approximate confidence interval procedures were 


evaluated by using computer simulations which are described in the next 
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chapter. During this evaluation process, the degrees of freedom in the 
expression ¥2>;« in equation (3.17), was increased and decreased from the 
defined values of F* given by this equation. The purpose of these mod- 
ifications was to find more accurate lower confidence limit procedures. 
The specific increases and decreases are described in chapter IV. The re- 
sults show that for some cases the procedure with modified degrees of 


freedom is more accurate. 


B. PARALLEL SYSTEMS 

An active parallel one-out-of-k system is defined as a system consisting 
of k subsystems such that system failure occurs when and only when all 
k subsystems fail. Equivalently, the system is successful when at least one 
of its subsystems is successful. Such a parallel system is said to be an ac- 
tive redundant system of order k. It can be shown that all of the prop- 
erties given in the preceding section for serial systems can be dualized to 
the corresponding properties for parallel systems by replacing any event 
by its complementary event. 

Let the parallel system be made up of k independently operating 
components, each with an exponentially distributed failure time and fail- 
ure rate A; ,1=1,2,...,.k. Then the system reliability, R(t), can be written 


as afunction of 4;,  i=1,2,...,.k as follows : 
k 
eo) — | 1a [ 1 —exp(- jt; ) | (3.20) 
(= 


Using equations that have been derived in the preceding section; 1.e., 
in the equations (3.5),(3.6),(3.14),(3.15),(3.16) and (3.17), the correspond- 
ing equations for the approximate 100(1 — «) % lower confidence limit 


for the reliability of a parallel system is 


I5 


A k A " 
R(t) 7 («) =|- L [1 = exp( = Am Ulayiti ) | (3.218 


for exponential failure times. 
For Weibull failure times the approximate 100(1 — a) % lower confi- 


dence limit for the reliability of a parallel system 1s 
a k NX A ie 
R(Oy(a) =! - He EES exp ea ae (3.22) 


C. SERIES-PARALLEL SYSTEMS 

| Series-parallel or parallel-series systems come in many varieties. One 
example of an active series-parallel system 1s defined as a system that is 
comprised of k subsystems in series connected to d subsystems in parallel 
or vice-versa. The reliability block diagram of an example system can be 


seen in Figure 3.1. 





Fig.3.1 Block Diagram Of Series-parallel System 


The system reliability R(t) of this series-parallel system can be written 
R(t) = L exp(— - At) JL1— TL {1 — exp(- 4,t)]] (3.23) 

if the failure time of each component Is exponential, and 
R(t) = exp(— z ePIC — T(t — exp(— 274P I] B.24) 


if the failure time of each component has a Weibull distribution. 

Using the same approach as those used for series systems and parallel 
systems for estimating the parameters, the corresponding equations for 
an approximate 100(1—«) % lower confidence limits for the reliability 


of the series-parallel system are 


A A r A 
R(t) 1(0) = L exp( — Aim, U(a) 2 it i (3.25) 


4 A m 
(=<T Ul —exp(—Ap ogy Pt 


_— Py 
and 


A _ 
Rr = L exp — An,u(a) 2, Peet yy, (3.26) 


A 
x 


(1ST [1 —exp(— 4m,u) Fal I 


IV. COMPUTER SIMULATIONS 


TEST PLAN : TESTING n, COMPONENTS UNTIL f, FAIL (RETP) 

RETP is a program written in FORTRAN, on the Amdahl mainframe 
computer. It performs the computer simulations that assess the accuracy 
of the lower confidence limit procedure for system reliability. A doc- 
umentation of this program and its associated subroutines is included in - 
Appendix C. 

The program accepts input parameters via an input file INI.DAT. 
For each replication, it generates the failure times for all the component 
items included in the test plan using a uniform random number generat- 
ing subroutine LRNDPC. The program determines the total test time 
accumulated for each component in the system and computes the esti- 
mates of the key parameters and the consequent lower confidence limit 
for system reliability for that replication. This is done for various system 
configurations; 1.e., series systems (Fig. 4.1), parallel systems (Fig. 4.2), 
parallel-series systems (Fig. 4.3) and a more general series-parallel system 
(Fig. 4.4). For each specific system configuration and set of input pa- 
rameters, the process is repeated 1000 times. When all replications are 
done, the routine EVAL processes the lower confidence limit estimates 
from all 1000 replications and determines the two measures of accuracy 
for the run, namely RSLOW and LEVEL. 

RSLOW is the 100(1 — «) percentile of the ordered set of lower con- 
fidence limits from the 1000 replications computed in a run. The true 
reliability of the system is RS. The closer RSLOW is to RS, the greater 


the accuracy of the procedure under evaluation in the run. If the proce- 


dure is exact, RSLOW will be equivalent to RS. To be conservative, 
RSLOW should be lower than RS. 

LEVEL measures the proportion of 1000 lower confidence limits, 
from a run with 1000 replications, which are /ower than the true system 
reliability RS. The closer LEVEL is to the specified confidence level for 
the procedure, | — a, the better the procedure. Values of LEVEL greater 
than (1 — a) reflect an under-estimation of RS which is conservative. 
Values of LEVEL less than | — & signal an over-estimation of RS which 


may be undesirable. 


W-Y)-OOO-OO-D)— 


Fig. 4.1 Block Diagram Of Series System 





Fig. 4.3 Block Diagram Of Parallel-series System 
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Fig. 4.4 Block Diagram Of General Series-parallel System 


(Rs) — 


wi 
us 


Simulation runs are performed using RETP for all combinations of 
failure time distributions and levels of key input parameters listed below. 
1. System 
a. 8 Exponential (Exp) components in Series (Case 1) 
b. 8 Weibull (Wei) components in Series (Case 2) 


c. 4 Exp and 4 Wei (Mixed) components in Series (Case 3) 
Results from these simulation runs are put into tabular form. 
By observing the most accurate result for different degrees of 
freedom, an approximate equation for obtaining the Degrees of 
Freedom can be derived, see Appendix B. Using this equation we 
then performed simulation runs for other systems for RETP. 


d. 8 Exponential components in Parallel (Case 4) 

e. 8 Weibull components in Parallel (Case 5) 

f. 4 Exp and 4 Wei (Mixed) components in Parallel (Case 6) 
g. 8 Exponential components in Parallel-series (Case 7) 
h 


. 8 Weibull components in Parallel-series (Case 8) 


ps | 


i. 4 Exp and 4 Wei (Mixed) components in Series-parallel (Case 9) 

j. 5 Exp and 10 Wei (Mixed) components in Series-parallel (Case 10) 
2. True System Reliability (RS) 

a. Hi (greater than 0.9) (Type A) 

b. Lo (greater than 0.8) (Type B) 


See pages 9-11 for description of the values of 1; and f,, for each of 
the above 10 cases. 


3. Level of Significance ( « ) 
a. 0.1 
b. 0.2 


4. Degrees of Freedom (DF) for the x? statistic as a function of the total 
number of failed test components (NFC) and total number of system 
components (NCOMP), (For cases 1,2 and 3 only). 


eh (Dy IN LC 


b. DF = 2(NFC + NCOMP) 
c. DF = 2(NFC - NCOMP) 
d. DF = 2NFC - NCOMP 


For case 4 up to 10 the equation derived in Appendix B was 
used. This equation is 


e. DF = 2NFC + 0.5(NC/NF) 
5. Test Plan for each component. 
a. Test 10 until 10 failures 
b. Test 15 until 15 failures 
. Test 15 until 11 failures 
. Test 15 until 7 failures 


a © 


e. Test 15 until 3 failures 


For the 8 exponential components in series in Case 1, the mission time 
for each of the components was chosen to be 10 hours. The chosen values 


of the scale parameters, 7, , are different depending on whether the system 
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is highly reliable (Type A system) or one with a lower reliability (Type B 
system). The ratios between the largest and the smallest failure rate were 
chosen to be 8 and 4.5 respectively for Type A and Type B systems. 

For the 8 Weibull components in series in Case 2, the mission time 
for each of the components was chosen to be 10 hours. In order to obtain 
a highly reliable system (Type A) or lower reliability system (Type B) the 
scale parameters, 4,, were chosen differently. The ratios between the larg- 
est and the smallest failure rate were chosen to be 8 for both system types. 
The shape parameter, f,, was chosen to be 2 for all subcases. The program 
will accomodate any value greater than zero for the shape parameter. 

In Case 3, which is a mixture of exponential and Weibull compo- 
nents, the mission time for each of the components was chosen to be 
equal to 10 hours. The scale parameter for each component type was 
chosen so that the ratios between the largest and the smallest failure rate 
were 4 for both component types. The shape parameter for each of the 
Weibull components is chosen to be 2 for all cases. 

In Case 4 ( 8 exponential components in parallel ), the mission time 
for each of the components was chosen to be 10 hours. The chosen values 
of the scale parameters, 4, , were different depending on whether the sys- 
tem is highly reliable (Type A system) or one with a lower reliability 
(Type B system). The ratios between the largest and the smalllest failure 
rate were chosen to be 1.525 and 2.575 respectively for Type A and Type 
B systems. 

For the 8 Weibull components in parallel in Case 5, no change was 
made for the mission time for each of the components as described in the 


previous paragraph. The ratios between the largest and the smallest failure 
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rate were chosen to be 1.28 for Type A systems and 1.63 for Type B sys- 
tems. The shape parameter, /,, was chosen to be 2 for all subcases. 

For a mixture of exponential and Weibull components in parallel, 
the mission time is still the same as before which is 10 hours. The values 
of the scale parameters, /;, for the exponential components were chosen 
to be equal to those for the Weibull components. The ratios between the 
largest and the smallest failure rates were chosen to be 1.3 and 2.2 for 
Type A and Type B systems respectively. The shape parameter, £,, for a _ 
Weibull component was chosen to be 2 for all subcases. 

In Case 7, the 8 exponential components in parallel-series (Figure 
4.3), the mission time for each of the components was chosen to be 10 
hours. The chosen values of the scale parameters, 7, were different de- 
pending on whether the system is highly reliable (Type A system) or one 
with a lower reliability (Type B system). The ratios between the largest 
and the smallest failure rate were chosen to be 8 for both Type A and 
Type B systems. 

For the 8 Weibull components in parallel-series (Figure 4.3) in Case 
8, no change was made for the mission time for each of the components 
which is 10 hours. The ratios between the largest and the smallest failure 
rate were chosen to be 8 for both systems. The shape parameter, f,, was 
set at 2 for all subcases. 

Case 9 is a mixture of exponential and Weibull components in 
parallel-series (Figure 4.3). The mission time for all components is 10 
hours. The values of the scale parameters, /,, for exponential compo- 
nents were chosen to be equal to those for Weibull components. The 


ratios between the largest and the smallest failure rate were chosen to be 
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4 for both component types and both systems. The shape parameter, f,, 
was set at 2 for all subcases. 

Case 10 is a mixture of exponential and Weibull components in a 
more general series-parallel configuration (Figure 4.4). The mission time 
for each of the components was chosen randomly. Similarly, the scale 
parameter for each of the components and shape parameter for each of 
the Weibull components were chosen randomly. This was done for both 
systems. 

Each simulation run of 1000 replications results in an output file 
OUT.DAT. The raw output from all the RETP runs are summarized in 
tabular form and placed in Appendix E. Each Table from Table 1A to 
Table 3B corresponds to a specific run case and system type combination 
using various degrees of freedom. Each Table from Table 4A to Table 
1OB corresponds to a specific run case and system type combination using 


only one specific degree of freedom. 


Pes 


V. ANALYSIS OF SIMULATION RESULTS 


After simulation runs were completed for selected cases, the results 
were analyzed for comparisons with results that had been obtained by 
Yee [ Ref. 7]. 

For each case, simulations were run with different degrees of freedom 
for the chi-square percentile point. This was done to determine if a for- 
mula for the degrees of freedom could be developed that would yield a 
more accurate lower confidence limit procedure. 

TEST PLAN: TESTING n, COMPONENTS UNTIL f, FAIL (RETP) 

When all components of the system have exponentially distributed 
failure times, the lower confidence limit procedures in this thesis are 
nearly the same as those developed by Yee. Consequently analysis of the 
simulation results will be discussed primarily for cases when some com- 
ponents of the system have failure times that have a Weibull distribution. 

A comparison of the four values of RSLOW, as in Table IA for each 
of 5 sampling plans (denoted by S/N), reveals that the lower confidence 
limit procedure with degrees of freedom equal to 2NFC - NCOMP 1s the 
most accurate lower confidence limit procedure. In S/N 2, for example, 
the RSLOW value of 0.9298 (using 2NFC - NCOMP degrees of freedom) 
is the largest such value below the RS value of 0.9305. The value of 
RSLOW above RS are optimistic and not as desirable as values of 
RSLOW which are equi-distant below RS. 

Table 2A and Yee’s Table 2A provide a comparison of the simulation 
results for eight Weibull components in series using the same key pa- 


rameters and the same mission times for each component as used by Yee. 
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When testing 15 components until all fail, Yee’s procedure and the 
procedure in this thesis have nearly equal accuracy. But when testing 15 
items until 7 fail, the procedure written in this thesis gives a more accu- 
rate value of RSLOW than the one from Yee’s procedure. In this case, 
Yee’s procedure yields values of RSLOW that are all above RS. In both 
Yee’s procedure and the one used in this thesis, the degrees of freedom 
of 27.NFC - NCOMP) gave the most accurate results among the four 
choices of degrees of freedom. It can also be seen that for more truncated 
testing, such as testing 15 items until 3 fail, the value of RSLOW tends 
to be higher and they are all slightly above RS. 

Table 3A displays the results of Case 3 for a type A system. The shape 
parameter, [, for the failure times with Weibull distribution was set equal 
to 2. The value of RSLOW resulting from the procedure developed in 
this thesis is more accurate than that for Yee’s procedure for various de- 
grees of freedom. Inspection of Table 3A reveals that the procedure 
corresponding to degrees of freedom 2(/NFC + NCOMP) its reasonably 
accurate for all 5 simulation cases and for both 90% and 80% confidence 
levels. 

Since the accuracy of the procedure depends on the extent of the 
truncation in the testing, a method was developed for choosing the de- 
grees of freedom that includes a term with the ratio NC/NF, where NC 
is the number of components placed on test and NF is the number of 
failed components. This makes the procedure dependent on the amount 
of truncation in the testing. Using information from all the simulations 
and the degrees of freedom that gave the most accurate values of RSLOW 
for all cases simulated, a formula for the degrees of freedom, DF, was 


developed using least squares methods. 


ar 


The equation is DF = 2NFC + I1.1(NC/NF) (see Appendix B for the 
derivation of this equation). Applying the equation to obtain RSLOW 
and observing the results for all cases simulated, resulted in a small 
modification to the above equation to yield the final equation DF = 
2NFC + 0.5(NC/NF). This final equation was used for the remaining 
cases that were simulated in RETP when any components had failure 


times with Weibull distribution. 
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Hable 1A : 8 Exp in Series, RS = 0.9305 (H1) 
min 2 = 0.0002 f/hr, max 2 = 0.0016 f/hr, UT = 10 hrs. 


SiN Test Deg. of 
“ Plan Freedom RSLOW LEVEL 
Test 10 until INFC ON 702 0.9630 
10 failed (160) 0.9256 0.9220 
NFC = 80 2(NFC + 0.9196 0.9970 
NCOMP 
(176) ) 0.9188 0.9910 
2NFC- 0.9296 0.9209 
NCOMP 
(152) 0.929] 0.8430 
tree | Gin | 0.9329 0.8400 


NCOMP 
(144) ) 0.9325 0.7300 











2 Test 15 until INFC 0.9277 0.9550 
IS failed (240) 0.9273 0.9080 
NEC = 120 a(n FC : 0.9233 0.9899 
(256) 0.9228 0.9750 
2NFC- ary 0.9298 0.9159 
NCOMP 
(232) . 0.9296 0.8329 
2NFC- | O1 | 0.9321 0.8440 
COMP 
(224) 0.2 | 0.9318 0.7470 
3 Test 15 until 2NFC -_ 0.9268 0.9550 


11 failed 


(176) 0.9262 0.9159 
2(NFC + 0.9208 0.9960 


NCOMP 
| (192) 0.9200 0.9880 
2NFC- — 0.9298 0.9159 


NCOMP 
(168) 0.9293 0.8430 
2(NFC- 0.9328 0.8430 


NCOMP 
(160) ) 0.9324 0.7350 


NFC= 88 


20 


Table 1A : 8 Exp in Series, RS = 0.9305 (H1) (Cont...) 
min 4 = 0.0002 f/hr, max 2 = 0.0016 f/hr, UT = 10 hrs. 


| Test Deg. of 
S/N 
Plan Freedom RSLOW LEVEL 
Test 15 until INFC 0.9241 0.9700 


i een (112) 0.9229 0.9310 
NFC= 56 2(NFC + 0.9145 0.9980 


NCOMP 
0198) 0.9130 0.9940 
2NFC- 0.9289 0.9190 


NCOMP 
| (104) 0.9280 0.8530 
2(NFC- 0.9338 0.8350 

NCOMP) 
(96) 0.2 0.933] 0.7200 
Con | Oot 0.9859 


(48) 0.9129 0.9750 
2(NFC + 0.8907 1.0000 


NCOMP 
(64) ) 0.8877 1.0000 
INFC- 0.9268 | 0.9439 


NCOMP 
(40) 0.9260 0.8600 
2(NFC- ar 0.9394 0.7530 


NCOMP 
Ps 0.9394 0.6339 

















+ 





Test 15 until 
3 failed 










NFC= 24 
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Yee’s results 
iar lA Ss Expin semes, KS = 0.931 (Hi) 
min A = 0.0002 f/hr, max A = 0.0016 f/hr, UT = 10 hrs. 


Test Deg. of 
S/N 
Plan Freedom RSLOW LEVEL 
Test 5 until INFC 0.919 0.982 


> vasa (80) 0.919 0.960 
NFC= 40 A~NFC+ | 01 | 0.906 1.000 


NCOMP 
(96) ) 0.905 0.999 
2NFC- 0.927 0.949 
NCOMP 
. (72) 0927 0.880 
7ore | GN) 0.934 0.821 
NCOMP 
(64) 0.934 0.702 
INFC 0.928 0.955 


(240) 0.927 0.908 
2(NFC + 0.923 0.990 


NCOMP 
(256) ) 0.923 0.975 
2NFC- 0.930 0.916 


NCOMP 
(232) 0.930 0.833 
2(NFC- 0.932 0.844 
NCOMP) 
(224) 0.932 0.747 
3 Test 15 until INFC mor 0.927 0.955 


M failed | (176 09260916 
NFC=88 A~ANFC+ | 0.1 | 0.921 0.996 


NCOMP 
(192) ) 0.920 0.988 
aNFc- | oa | 0930 | 0916 
NCOMP 
(168) 0.929 0.843 
2ANFC- | 0.1 | 0.933 0.843 
aaa We 0.932 735 
(160) : 0. 








Test 15 until 
15 failed 






NFC= 120 








a 


Yee’s results 
Table 1A : 8 Exp in Series, RS = 0.931 (H1) (Cont...) 
min A = 0.0002 f/hr, max A = 0.0016 f/hr, UT = 10 hrs. 


Test Deg. of Measures of accuracy 
S/N 
Plan Freedom RSLOW LEVEL 
Test 15 until INEC | Fon 0.924 0.970 


{| ails (112) 0.923 0.931 
NECG=56 °NFC+ | 6.1 | ¥ olois 0.998 


NCOMP 
COMP) ee a 
INFC om fos 


NCOMP 
NEG? | OD Wiss4 0.835 


a 
Test 15 until NFC | -O | 915. | aoe 
NFC=24 AxNFC+ | ot | 7 0.891 |__G60nay 


NCOMP 
(64) ) 0.888 1.000 
ee, | Om. Wo 0.944 


NCOMP 
(40) 0.926 0.860 


2(NFC- 0.939 0.753 
NCOMP 
(32) ) 0.939 0.634 
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Table 2A : 8 Wei in Series, RS = 0.9798 (H1) 
min 2 = 0.001 f/hr, max 7 = 0.008 f/hr, UT = 10 hrs. 


Test Deg. of 
S/N 
Plan Freedom RSLOW LEVEL 
Test 10 until ONFC | 0.1 | 0.9602 0.9859 


10 failed (160) 0.9506 0.9820 
NEC = 80 2(NFC + 0.9565 0.9890 


NCOMP 
“tie” 0.9461 0.9870 
ae | Oy osers 0.9820 
NCOMP 

(152) 0.9529 0.9809 
2ANFC- | 0.1 | 0.9638 0.9800 
NCOMP) 

(144) 0.9553 0.9770 
Test 15 until INFEC 0.1 0.9706 0.9770 


esse (240) 0.9649 0.9719 
NFC= 120 ANFC+ | 0.1 | 0.9687 0.9820 


NCOMP 

(256) 0.9627 0.9790 
2NFC- | 0.1 {| 0.9715 0.9730 

NCOMP 
(232) bee 0.9659 0.9640 
2(NFC- 0.9724 0.9650 

NCOMP) 
(224) a 0.9671 0.9590 
Test 15 until INFC [isle 0.9704 0.9719 


it ies (176) 0.9637 0.9650 
NFC=88 2NFC+ | 0.1 | 0.9679 0.9800 


NCOMP 
(192) > 0.9606 0.9760 
2NFC- 0.9716 0.9660 


NCOMP 
(168) 0.9652 0.9590 
ANFC- | 0.1 | 0.9729 0.9590 


NCOMP 
oN 0.9668 0.9510 
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Table 2A : 8 Wei in Series, RS = 0.9798 att (Cont...) 
min A = 0.001 f/hr, max 2. = 0.008 f/hr, UT = 10 hrs. 


[wT oe | es | 
ae Freedom 
Test 15 until DAC 

— ie) 

NFC= 56 2NFC + 

2NFC- | 0.1 | 0.9779 | 0.9240 


NCOMP 
(104) 0.9691 0.9110 
2(NFC- 0.9795 0.9080 


NCOMP 
(96) 0597 53 0.8900 
INFC 0.9889 0.8000 
(48) 0.9814 0.7770 
2(NFC + 0.9857 0.8469 


NCOMP 
(64) ) 0.9758 0.8360 
2NFC- 0.9906 0.7510 


NCOMP 
(40) 0.9843 0.7280 
2(NFC- 0.9922 0.6890 


Re 
(32) ) 0.9872 0.6530 
































Test 15 until 
3 failed 











NFC= 24 
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Yee’s results 
Table 2A : 8 Wei in Series, RS = 0.980 (Hi) 
min A = 0.001 f/hr, max 4 = 0.008 f/hr, UT = 10 hrs. 


Test Deg. of 
Plan Freedom RSLOW LEVEL 
Test 5 until INFEC 0.947 Oe 


SHANES (80) 0.930 0.989 
NFC = 40 2NFC+ 0.937 0.994 


NCOMP 
(96) 0.918 0.993 
2NFC- 0.951 0.989 

NCOMP 
(72) 0.937 0.986 
2(NFC- 0.956 0.985 

NCOMP 
(64) ) 0.943 0.981 
INFC Pacaie- | 0.978 0.918 
(240) 0.974 0.913 

ANFC+ | 0.1 | 0.977 0.931 

NCOMP 

. (256) ) 0.972 0.924 
2NFC- o979 =| tA 

NCOMP 

(939) 0.975 0.901 
2(NFC- 0.980 0.904 

NCOMP 
(224) 0.975 0.889 

3 Test 15 until INFC 0.982 0.876 

il failed (176) 0.977 0.860 

NFC= 88 2NFC : 0.980 0.894 

(192) 0.975 0.882 

2NFC- 0.983 0.861 
NCOMP 

(168) 0.978 0.839 

2Nre: yor 0.983 0.840 
NCOMP 

(160) 0.979 0.819 








Test 15 until 
15 failed 






NFC= 120 
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Yee’s results 
Table 2A : 8 Wei in Series, RS = 0.980 (H1) (Cont...) 
min 2 = 0.001 f/hr, max 2 = 0.008 f/hr, UT = 10 hrs. 


S/N Test Deg. of 
[i Plan Freedom RSLOW LEVEL 
Test 15 until NEG. WO ee 0.800 


— (112) 0.981 0.779 
NFC = 56 A~NFC+ | 0.1 | 0.985 0.839 


NCOMP 
(128) 0.978 0.824 
2NFC- 0.988 0.776 
NCOMP 
(104) 0.982 0. a8 
2(NFC- 0.989 0.746 
NCOMP 
(96) ) a 0.983 O32 
Test 15 until ONFC -- 0.1. | 0.994 0.621 


9 tailed (48) 0.991 0.584 
NFC=24 ANFC+ | 0.1 | 0.993 0.705 


NCOMP 
(64) ) 0.988 0.685 
2NFC- | 0.1 | 0.995 0,548 


NCOMP 
(40) 0.992 0.514 
2°NFC- | 0.1 | 0.996 0.408 


NCOMP 
(32) ) 0.993 0.417 
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Table 3A : 4 Exp and 4 Wei (mixed) in Series, RS = 0.9801 (H1) 
min A = 0.0002 f/hr, max 1 = 0.0008 f/hr, UT = 10 hrs. 


SIN Test Deg. of 
y Plan Freedom RSLOW LEVEL 
Test 10 until rea, 0.1ib.. 90.9796 0.9190 


Jeter (160) 0.9791 0.8720 
NFC= 80 2(NFC + 0.9777 0.9740 


NCOMP 

(176) 0.9771 0.9489 

2NFc- | 0.1 | 0.9805 0.8720 
NCOMP 

(152) 0.9801 0.7980 

2(NFC- 0.9815 0.7950 
NCOMP 

is DE 

2 Test 15 until in OU Te 0.8950 


1S failed (240) 0.9799 0.8250 
NFC= 120 ANFC+ | 0.1 | 0.9789 0.9510 


NCOMP 
(256) ) 0.9786 0.9209 
2NFC- | 0.1 | 0.9808 0.8370 

NCOMP 
(232) 0.9805 0.7530 
ANFC- | 0.1 | 0.9814 0.7670 

NCOMP) 
(224) 0.9812 0.6670 
3 Test 15 until INFC | 0.1 | 0.9802 0.8920 


Lp ules (176) 0.9797 0.8400 
NFC=88 2(NFC+ Oi 0.9785 0.9590 


NCOMP 
(192) ) 0.9779 0.9320 
2NFC- | 0.1 | 0.9810 0.8410 


NCOMP 
(168) 0.9806 0.7550 
2ANFC- | 0.1 | 0.9819 0.7580 


NCOMP 
ve 0.9815 0.6590 
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Table 3A : 4 Exp and 4 Wei in Series, RS = 0.9801 (Hi) (Cont...) 
min 2 = 0.0602 f/hr, max 2 = 0.0008 f/hr, UT = 10 hrs. 


SIN Test Deg. of 
is Plan Freedom RSLOW LEVEL 
Test 15 until INFC | 0.1 | 0.9803 0.8950 


Hi feulles! (112) 0.9797 0.8260 
NEC = 56 2NFC+ | 0.1 | 0.9777 0.9690 


NCOMP 
(128) ) 0.9770 0.9439 
INES... | mole Mewtiosis 0.8060 


NCOMP 
(104) 0.9811 0.7300 
ANFC- | 01 | 0.9828 0.7090 


NCOMP 

(96) ) 0.9825 0.6000 
ome | On Goss 0.8730 
(48) 0.9800 0.8020 

ANFC+ | 0.1 | 0.9757 0.9740 

NCOMP 

64) ) 0.9739 0.9560 


2NFC- 0.9839 0.7250 


NCOMP 
(40) 0.9831 0.6260 
ANFC- | 0.1 | 0.9868 0.5150 


NCOMP 
32) ) 0.9862 0.4140 




























Test 15 until 
3 failed 





NFC= 24 
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ees Tesults 
Table 3A : 4 Exp and 4 Wei (mixed) in Series, RS = 0.980 (H1) 
min 7 = 0.002 f/hr, max 7 = 0.008 f/hr, UT = 10 hrs. 


Pv Test Deg. of 
i Plan Freedom RSLOW LEMEL 
Test 5 until wFc | 01 | 0.979 0.942 


> failed (80) 0.978 0.905 
NEC=40 ANFC+ [{ O01 | 0.975 0.987 


NCOMP 
(96) 0.974 0.976 
res | 01 0.981 0.881 
NCOMP 
(72) 0.980 0.805 
2(NFC- 0.983 0.771 
NCOMP 
(64) L 0.982 0.684 
2 Test 15 until NEG | COl  Oasl 0.863 
1S failed (240) 0.980 0.800 
NEC = 120 2NFC F | Cn. “O75 0.94] 
86)» 0.979 0.898 
2NFC. | Ol | 0.981 | 0.881 
NCOMP 
(939) 0.980 0.805 
ANFC- | 01 | 0.982 0.725 
NCOMP 
(224) ) 0.981 0.631 
2NFC |_O.1 | 0.981 0.864 


a) 0.980 | 0.801 
ANFC+ | 0.1 | 0.979 0.951 


NCOMP) 
(192) 0.978 0.907 
2NFC- | 0.1 | 0.982 0.802 
NCOMP 
(168) 0.981 0.698 
ANFC- | 0.1 | 0.982 0.702 
NCOMP 
(160) 0.982 0.591 











Test 15 until 
11 failed 






NFC= 88 


3° 


Yee’s results 
Table 3A : 4 Exp and 4 Wei in Series, RS = 0.980 (Hi) (Cont...) 
min 7 = 0.002 f/hr, max 4 = 0.008 f/hr, UT = 10 hrs. 


| Test Deg. of 
S/N 
Plan Freedom RSLOW LEVEL 
Test 15 until INFC | 0.1 | 0.981 0.865 


7 failed (112) 0.980 0.787 
NFC= 56 2(NFC + eee Cee 


NCOMP 
128) 0.978 0.920 
2NFC- | 0.1 | 0.982 0.769 


NCOMP 
(104) 0.982 0.676 
71GN1 cL On CO El Osx 0.644 


NCOMP 
(96) 0.983 0.523 
Test 15 until IFC | 01 | 0.982 0.843 
ape (48) | 02 | 0.981 | _—0.762 
NFC=24 ANFC+ | 0.1 [| 0.976 0.970 
Nee? aoe 0.975 
(64) : . 


Nec. | OuIN Woow 0.684 


NCOMP 
(40) 0.984 0.580 


2(NFC- 0.987 0.459 
NCOMP 
32) 0.987 0.386 
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Table 4A : 8 Exponential in Parallel, RS = 0.9345 (Hi) 
min 2 = 0.1000 f/hr, max 2 = 0.1525 f/hr, UT = 10 hrs. 
Degrees of Freedom = 2NFC 


yN | Test wip. Of 
— Freedom RSLOW LEVEL 
Test 10 until | 0.1 | 0.9306 0.9290 


10 failed 
160 
0.9320 0.8180 


Test 15 until veel ane 0.5 0.9190 


15 failed 740 















NFC = 80 


NFC= 120 


0.2 0.9324 0.8210 


F Test 15 until | 0.1 | 0.9309 0.9240 





11 failed 176 
0.2 0.9312 0.8329 
NFEC= 88 


0.9307 0.9209 


112 
37 0.9308 0.8310 
NEG = 56 


Test 15 until 0.9280 0.9220 


3 failed 
0.2 0.9286 0.8370 
NFC=24 


Test 15 until 
7 failed 







4] 


Table 5A : 8 Wei in Parallel , RS = 0.9265 (Hi) 
min A = 0.100 f/hr, max 2 = 0.128 f/hr, UT = 10 hrs. 
Degrees of Freedom = 2NFC + 0.5(NC/NF) 


Test ee, et 
Plan Freedom RSLOW LEVEE 


Test 10 until | 01 | 0.9154 0.9540 


10 failed 
0.2 Veo, 0.8880 
NFC= 80 


Test 15 until | 0.1 | 0.9179 0.9510 


15 failed 
0.2 0.9170 0.8770 
NFC= 120 


Test 15 until 0.9248 0.9180 


11 failed 
0.2 0.9250 0.8130 
NFC= 88 


Test 15 until 0.1 | 0.9323 0.8630 


7 failed 
0.2 0.9318 0.7589 
NFC= 56 


Test 15 until po oe 0.9270 


3 failed 
0.2 0.8980 0.8850 
NFC= 24 
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Table 6A : 4 EXP and 4 Wei (mixed) in Parallel, RS = 0.9408 (Hi) 


min A = 0.100 f/hr, max 4. = 0.130 f/hr, UT = 10 hrs. 
Degrees of Freedom = 2NFC + 0.5(NC/NF) 


S/N Test Deg. of 
Plan Freedom RSLOW LEVEL 


Test 10 until 
10 failed 


NFC= 80 


Test 15 until 
15 failed 


NEC = 120 
Test 15 until 
11 failed 
NFC= 88 
Test 15 until 
7 failed 
NFC= 56 


Test 15 until 
3 failed 





NFC= 24 











Measures of accuracy 










Prony Woisasah 2019240 
Vor 0.946 |__ 0.8920 







rol ip ~ 0383 0.9080 


SI] 
0:93:25 0.8410 
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Table 7A : 8 Exp in Series- Parallel, RS = 0.9249 (Hi) 
min A = 0.0003 f/hr, max 4 = 0.0024 f/hr, UT = 10 hrs. 
Degrees of Freedom = 2NFC 


Test Deg. of <a of accuracy 
le Freedom / RSLOW | LEVEL 


Test 10 until | 0.1 | 0.9241 0.9159 


10 failed 
_ 0.9226 0.8620 


Test 15 until Ra os 0.8960 


15 failed 
0.2 0.9236 0.8360 
NFC= 120 


Test 15 until | 0.1 | 0.9252 0.8950 


11 failed 
0.2 0.9226 0.8390 
NFC= 88 


Test 15 until or | woo 0.9260 


7 failed 
0,2 0.9204 0.8789 
NFC= 56 


Test 15 until | OL | 0.9153 0.9620 


3 failed 
; 0.9123 0.9260 





44 


Table 8A : 8 Wel in Series- Parallel, RS = 0.9328 (Hi) 
min A = 0.002 f/hr, max 4 = 0.016 f/hr, UT = 10 hrs. 
Degrees of Freedom = 2NFC + 0.5(NC/NF) 


Test Wee of 
Plan Freedom RSLOW LEVEL 


Test 10 until | 0.1 | 0.9106 0.9579 



















10 failed Al 
0.2 0.8939 0.9510 
NFC= 80 
2 Test 15 until Per | ooo 0.9330 





15 failed 





241 
0.9108 0.9230 


NFC= 120 





Test 15 until | 0.1 | 0.9266 0.9220 


11 failed 






0.2 0.9149 0.9119 
NFC= 88 


Test 15 until | 0.1 | 0.9414 0.8609 


7 failed 
0.2 0.9250 0.8419 
NFC= 56 


Test 15 until 0.9699 0.7079 


3 failed 
0.2 0.9536 0.6740 
NFC= 24 
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Table 9A : 4 EXP and 4 Wei in Series- Parallel, RS = 0.9276 (Hi) 
min 4 = 0.005 f/hr, max 2 = 0.020 f/hr, UT = 10 hrs. 
Degrees of Freedom = 2NFC + 0.5(NC/NF) 


S/N Test Deg. of Measures of accuracy 
Pat Freedom RSLOW LEVEL 


Test 10 until | 0.1 | 0.9096 0.9400 


10 failed 
0.2 0.8904 0.9349 
NFC= 80 


















Test 15 until Ol | Womess 0.9170 
15 failed 74) 
0.2 0.9079 0.9040 
NFC = 120 
Test 15 until | 0.1 | 0.9266 0.9050 
11 failed aE 
0.2 0.9109 0.8850 
NFC = 88 


Test 15 until or | Weleses 0.8550 


7 failed 
OF 0.9203 0.8390 
NFC= 56 


Test 15 until Poi | soseen 0.7430 


5 
3 failed 
0.2 0.9451 0.7200 
NFC= 24 





46 


Table 10A :10 EXP and 5 Wei in Series- Parallel, RS = 0.9472 (H1) 
Exp: min 2 = 0.025 f/hr, max 4 = 0.075 f/hr. 
Wei: mind = 0.055 f/hr, max 4 = 0.095 f/hr. 
UT = 10 hrs. 
Degrees of Freedom = 2NFC + 0.5(NC/NF) 


S/N Test DOE of 
Plan Freedom RSLOW Aue 
Test 10 until ee ORT 0.9710 
301 
peers a, 


10 failed 
Test 20 until | 0.1 | 0.9449 0.9370 


2 
20 failed 601 
0.2 0.9432 0.8950 
NFC = 300 
Test 20 until 0.9444 0.9349 










NFC=150 










15 failed AS] 
0.2 0.9430 0.8979 


0.9435 0.9290 


301 
avais | 0.9000 


0.9344 0.9809 


le2 
OF2 0.9310 0.9520 





DEG 225 


Test 20 until 
10 failed 











NFC= 150 


Test 20 until 
5 failed 






NFC=75 
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VI. CONCLUSIONS 


The lower confidence limit procedures developed in this thesis are 
modifications and extensions of the lower confidence limit procedures 
written by Yee. Procedures were developed and evaluated in this thesis for 
more complex system structures than Yee examined. In addition a dif- 
ferent method for estimating the parameters, in the Weibull distribution 
was uSed here rather than the maximum likelihood procedure used by 
Naser 

The evaluations show that the approximate lower confidence limit 
procedures are reasonably accurate for all system structures examined if 
the degrees of freedom are chosen judiciously. Although the equations 
given here for choosing an appropriate value for the degrees of freedom 
are for the system simulated in this thesis, it would be prudent to run 
simulations for complex systems that differ substantially from those ex- 
amined in this thesis in order to determine an appropriate number for the 
degrees of freedom in the lower confidence limit equation. 

The degrees of freedom equation derived in Appendix B, DF = 
2NFC + 0.5(NC/NF), yielded lower confidence limit procedures with 
reasonable accuracy for test plans that were truncated on the number of 


failures. 
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VII. APPLICATION EXAMPLES 


Based on the procedures evaluated by the RETP runs, three config- 
urations of systems, a specific test plan and failure time data sets were 
constructed to illustrated the use of the procedures in providing a lower 
1001.1 — a) % confidence limit for system reliability. Actual results of one 
computer run for each case simulated are given below. 


8 Exponential components in Series 
TEST PLAN 1 - Test 15 until 7 fails for each component 


I. Raw Data 


Com Ordered Failure Times 
al T(1) T(2) T(3) T(4) cS ) T(6) a7 ) 
1 695.0 Zl 3 1965.7 2628. 8 3304. 9 3946. 3 4014. 2 
2 e255 Dee | 230. 8 Bo 2 401.4 1071.9 1222-0 
3 ii. 9 394.7 G22 6 506.5 519.5 520. > 582.6 
4 525..5 B50. 2 Fa) Fi 837.8 844.2 873.7 894.4 
5 62. 6 110.0 12423 126.8 325. & 384. 0 502.6 
6 20. 0 85.6 107.0 108.0 161.5 190.1 201. 6 
7 ae), 77.9 100.4 ee. 156.0 180.9 193.5 
8 Ge 1 34.6 48.1 65.3 95.2 136. 6 160.4 

II. Data Summary 

Comp(i) UTCi) NCCi) NFCi) AEC 1) ELM(i) ERC i) ERC i)*TTCi) 
ai 1070 15 y 49310. 2 0.00012 O° 037i10 1829.9 
2 10.0 15 7 13572.8 0.00045 0. 13683 1829.9 
3 10.0 ls fi Ti 57a 0.00077 0.23589 1829.9 
4 O20 15 7 11728. 8 0.00051 0.15601 1829.9 
5 10.0 5 i 5658.4 0.00106 0. 32339 1829.9 
6 10. 0 iS Vl 2487.0 0.00241 0. 73578 1829.9 
7 10.0 15 yi 2433. 2 0.00247 @275205 1829.9 
8 10.0 15 7 1829.9 0.00328 1. 00000 1829.9 


49 


RS 0.93054 


ALPHA Ore 
NFC 56 

CHISQD 104 122.86 (from the table Chi-square distribution ) 
LMU 0.00419 





RSLOW 0.92894 








CASE 4 : 8 Exponential components in parallel 


cserece TEST PLAN 1 - Test 15 until 7 fails for each component 

I. Raw Data 

Comp Ordered Failure Times 
a T(1) TC) TCS) T(4) T(5) T(6) TOD 
1 1. 3899 2.4825 2.7315 5.2576 46.6058 Fe (se) o3 7) 8.0285 
2 0.530355 0. 79a 0. 8586 127 1.4936 3.9886 4.5471 
3 0.5839 2.0591 2.2047 2.6425 2.7106 2.9156 GaGeee 
4 2.1258 2.3264 2.8844 5.4715 5.5130 5.70555 “Sasai@ 
5 0.4813 0.8465 0.9563... 019751) “2esg6s 2.9542 3.8674 
6 0.1748 0.7471 0.9340 £0.9424 1. 4090 Jeno 9 3 1.7599 
7 0.3296 0.7519 0.9697 lV 3716 1. 5065 1. 7469 1. 8685 
8 0.0638 0.3635 0.5050 £0. 6849 0.9989 1.4334 1.6832 

II. Data Summary 

Comp(i) UTCi) NCCi) NFCi) Tita ELM( i) ER(i) ERC) iG 
iL 10.0 15 7 98.6204 0.06084 0. 19467 19, 1986 
2 10.0 15 7 49.7592 0.12058 O2g3583 19.1986 
3 10.0 15 7 40.4724 0.14825 0.47436 T9M1986 
+ 10.0 15 yi 76. 5961 905.07 Ges 0.25065 19. 1986 
5 10.0 15 7 435263 “O.le785 0.44108 19, 1986 
6 TOE0 15 7 21.7044 0.27644 0.88455 19. 1986 
y, 10.0 13 Z 23.4926 0.25540 0.811 7 22 19. 1986 
8 10.0 15 y/ 19.1986 (Oveav52 1. 00000 19, 1986 


Parameter df Value 
RS 0.93459 
ALPHA 0.1 
NFC 56 

CHISQD po 131.56 (from the table Chi-square distribution ) 
LMU 0.42828 
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RSLOW 0. 93069 


Case 7 : 8 Exponential components in Series-Parallel. 
eeeee-- (4 components in parallel connected to 4 components in series ) 
cocece- TEST PLAN 1 - Test 15 until 7 fails for each component 


I. Raw Data 


Comp Ordered Failure Times 
i 1) TiG2) T(3) T(4) T(5) T(6) TG?) 
if fosweO9 627me05 9l0N497 1752.519 2203. 295 26304893 2676. 161 
2 Oooo Weeeteet? 1535637 210.102 267.609 714.622 814.687 
3 747604 2635110 281.713 337.648 346.349 372.313 388.416 
4 CivenOOs 2Zogea9) 2947454 556.552 562.782 582.436 596.273 
5) 41.716 73.360 624602 Semel) 2172230 256.028 335.172 
6 3. 354 27 a0e3 71.346 71.990 107.635 126.749 134.425 
7 TA aM se. aie Sy 66.954 94/03 V104e022 120.616 129.016 
8 4.053 Zo 8095 322090 43.517 63.475 oie oO (106.955 


Comp(i) UTCi) NCCi) NFCi) aC 1) ELM(i) ER(i) ER(i)*TT(i) 
1 10.0 ILS: 7 32873.465 0.00018 0.03711 E219. 9089 
Zz TOE 0 ILS. 7 8915.184 0.00067 0.13683 a9. 9089 
3 10.0 15 7 Si7 27 7 0 0G116.. 0.23589 i219..9069 
4 10.0 Ls 7 T8197 Sie 0700077 0. 15601 ag Ge? 
5 10.0 i 7 S7e gee” O. U0ISIR 07525359 E219 2089 
6 102.0 Ie 7 1657. 97a 0.00362 0.73576 1219. 9069 
7 1O%0 ES 7 io227 109° 0. 00870 “0275205 W219. 9069 
8 10. 0 15 7 719. 909 7 0.00492 i. 00000 1219. 90e0 


Parameter df Value 
RS 0. 92496 
ALPHA 0.1 
NFC 56 
CHISQD a2 131.56 (from the table Chi-square distribution ) 


LMU 0.00674 











RSLOW On2223 


Case 2: 8 Weibull components in Series 
-—<eee-~ TEST PLAN 1 - Test 15 until 7 fails for each component 
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Raw Data 


Ordered Failure Times 
T(3) 


T(2) 


TC 


4) 


A Rs) Joa Ob 
. 43E+05 
5h og eal 
69E+10 
29n+-05 
25 tev 
a BO 7 
.41E+04 


T(5) T(6) T(7) 
406.506 444.204 448.009 
100.177 163.702 174.788 

93.052 96.477 98.541 
102.724 104.502 105.736 

57.083 61.971 70.905 

36.680 39.804 40.992 

33.384 35.949 37.179 

24.394 29.221 31.666 
i) ER(i)  ER(i)*TT(i) 
-07 0.23E-04  4107.8930 
-03 0.96E-01 4107.8930 
-09 0.33E-06 4107. 8930 
-08 0.60E-06  4107.8930 
-03 0.14E#00 4107. 8930 
-04 0.17E-01 4107. 8930 
-05  0.27E-02 4107. 8930 
-02 0.10E+01 4107. 8930 


(from the 


table Chi-square distribution ) 


Ordered Failure Times (h) raised to the power of Beta 
DGG) 


Tt (ip T 2) 
.14E+07 0. 30E+07 
.68E+03 0. 94E+03 
25E+08 0. 42E+09 
64E+08 0.77E+08 
.37E+03 0. 62E+03 
.94E+03 0. 66E+04 
.77E+04 0. 30E+05 


Gs) 


ooo 0 00 © 


SSD E+07 
. LOE+04 
. 49E+09 
~ IZEs09 
Rao} 12 1015, 
. 89E+04 
» 45E+05 


T' (4) 


oOoo0o0o0 Oo 


a2 


. 84E+07 
. 12E noe 
. 74E+09 
. 49E+09 
- 7 PEs 
. 9OE+04 
. (9Eaw> 


oR eee kexeke’ 


T' (5) 


. IEEE OS 
. ISEFOS 
S709 
Re gOS) 
. 17E+04 
el SD ol OB 
Ge) lagcUi5) 


ooo 0d O00 


. 14E+08 
. 34E+04 
. JSEROO 


54E+09 


. 20E+04 
. LOBES 
- E06 


CO Oe Oo oO cee 


. 15E+08 
. 38E+04 
. LOEREO 
. > EO? 
» 25E+04 
. ZTErs 
. 1SEe6 


8 Ope O20. OGEte2eee 1 SET+OS) 0. 16E4O3 0.23Et03 O0.31Et03 0. 35E+03 


Case 5: 


8 Weibull components in parallel. 
TEST PLAN 1 - Test 15 until 7 fails for each component 


I. Raw Data 


Comp Ordered Failure Times 

a m1) TG2) T@2) T(4) Gey) Ge) T(7) 

iL. 3.726 4.982 yg 4.¢ Peco 8.130 8. 884 8. 960 

2 2. 296 2,005 z. W2¢ 3.414 3.48 55 6.296 Ge 7235 

3) Z. 390 4.506 4.662 5.104 5.070 5. 360 5.475 

4 4,556 4.766 3. 307 72510 7.837 7.464 (Seas) 

5 Ze i156 2. 860 3. 040 3. 069 4.921 5. 342 Ow lio 

6 ibe Zeo7 1. 2.986 3. 000 3.668 3. 980 4.099 

7 1. 763 2. 663 3.024 3.520 5 9 4,059 4.198 

8 0.771 le659 2.168 Zo 3.049 32655 Jno56 

II. Data Summary 
Comp(i) UT(i) NCCi) NF(Ci) TCL) ELM(i) BRGi) WREH)*TI( i) 

i OBES, is) 7 0. 45E+04 Omler=02 0e27E-0)) 120.5522 
Z 10.0 15 7 0. 24E+03 ONSOEe0l SOMSdE+00) 120.5522 
3 LOZ 1S 7 0. 2 705 OnZ26GRe0S mOw45E=02 120.5522 
4 10.0 i) 7 0. 76E+05 Uae oae Oe lobe02 120, 5522 
5 kOe 16 15 7 OS Se 01} QO2226b-01 0. 388700, 120.5522 
6 10-0 a) 7 0252E+03 Oeber-Ol 02238400. 120.5522 
7 1020 15 7 0. 13E+04 O-55be02 0.95E501 120.5522 
8 10.0 1S 7 0. 1ZE40S 0. 586-00) O.MMOE+OL 020.5522 





(from the table 


Parameter df Value 
RS 0.92659 
ALPHA Omi 
NFC 56 
CHISQD 116 135. 89 
LMU 0.07046 
RSLOW 0.92598 


Chi-square 


distribution ) 


Ordered Failure Times (h) raised to the power of Beta 
i [Clg nee) LaeSy, eece ) TO) i Co) T @ 


5) 


il 0.35E+02 0. 77ETO2 0. 88ETOZ™0. 21E+US" 0. 29E+05 07 7 ETO See Oo 
2 0. 38E+O1 0.52E+01 0.55E401 0. 7IETO] 0. 86ET0OT 0, MORaO NCR eer eZ 
3 0.52E+02 0.90E+03 0.10E+04 0.16E+04 0.17E+04 0. 20E+04 0. 22E+04 
4 O. 7OE+03 0. 86E+03 0. 14E+04 0.54E+04 0.55E+04 0. 60E+04 0. 63E+04 
5 O.41E+01 0. 69E+01 0. 77E¥FO1 0. 7EEFO1 UO. I9ETO2 0. 22h 702m Er o2 
6 O. 20E+O1 0. 14E+02 0. 19E4#0O2 0. 19E+02 OF 32E+02 ©. GOETO2 OFeaET0Z 
7 0.63E+01 0.24E+02 0.37E+02 0. 65E102 0: 7/5E1+02 0. 96E+029ORii ee 
8 0.64E+00 0. 28E+01 0.37E+01 0.48E+01 0.66E+01 0.90E+01 0. 10E+02 

Case 8 : 8 Weibull components in series-parallel. 

---==—-= (4 components in parallel connected to 4 components in series) 

-—- === TEST PLAN 1 - Test 15 until 7 fails for each component 

I. Raw Data 

Comp Ordered Failure Times 
i TOL) AC TG) T(4) EGS) T(6) 7) 
i 5.142 6.872 Bs IAYS, 10.001 1, 214 2. 254 12, 353 
Z 3.184 By 81802, 4.051 4.734 Deas 5. 731 O322 
33 3. 344 62279 6.497 7,118 7.204 7.469 7.022 

4 Omer 9 62.07 3 7 30 10. 234 10.272 10. 450 10.574 
5 83.383 110.574 117.532 18666] 190.2765 2064571, Zig 
6 25.840 53.419 59.727 59.986 —(735300° 979.6038  "oimeae 
7 24.289 36.688 41.663 49551 55175931 j55.920 Syme 
8 Gag i, 619 23.126 26.9319) 32.526 38.962 eee 

II. Data Summary 

Comp(i) UTCi) NC(i) NFCi) TEG® ELM(1) ERG?) ER(i)*TI(i) 
1 ILO 0 15 Z O.11E+05 0.66E-03 0.38E-0O1 39955 
Zz 1 08, 15 } 0.40E+03 0. 18E=0) O00 10E+0 S399%scer. 
3 10.0 iS 7 0.12E+06 0.58E-04  0.33E-02 3995665 
4 FOn.O 15 7 0.33E+O06 0O.2Z1E-06 0.12E-02 39973508 
5 10.0 1S 7 0.26E+06 O0.27E-04 0.15E-02 39905eu8 
6 10.0 als) 7 0. 16E+07 0.44E-05 0.25E-03 33805008 
i 10.0 rS 7: 0.65E+0O7 O.11E-05 0.61E-04 39935008 
8 10.0 i> i; 0.67E+O4 0.10E-02 0.60E-01 3990ereus 


Parameter folic Value 
RS 0.93682 
ALPHA Oi 
NFC 56 
CHISQD 104 122.86 (from the table Chi-square distribution ) 
LMU 


0.01922 


ee ei ee ee ee ee 
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RSLOW 


SS 


IV. Workarea 


0. 93382 


Ordered Failure Times (h) raised to the power of Beta 
T'(1) Ta@) 
. 84E+02 
Bose 1 


23E+03 
30E+04 
34E+04 


. 60E+04 
-ooErOS 
BoOLtU2 


Paey 


EOS 
Pov Fol 
. 40E+04 


37E+04 


. 57E+04 


42E+05 


| IS E06 
. 16E+03 


M21E+03 
Soon re | 
-47E+04 
. 58E+04 


63E+04 


OS 
SIE T06 
« 2 EOS 


T'(4) 


. ouiO 3 
. ZO 2 
. 71E+04 
2 UD 
- 65E+04 
. 98E+05 
O00 
27 id-OS 


B (5) 


. 69E+03 
. 1GE+02 
» 75E+04 


24E+05 
15E+05 


elie 1S) 
. 39E+06 
. 37E+03 


ied) 


26E+05 
18E+05 


. 12E06 
. 49E+06 
Jo OleOS 


e9OE+03 
eo OZ 


97E+04 
27 oO 
Z35E705 


5 ls de els: 
yeep sils 
oo OS 
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APPENDIX A. DERIVATION OF FORMULA FOR BIAS VALUE 


Suppose Y has an extreme-value distribution; 1.e., Y~ EVo, yp). 
Suppose Voss) S Yar, S Yous) SS Yyig are the middle n-s ordered 
Statistics in a random sample of size n from extreme-value distribution 


with pdf. and Cdf. given by 


Hosea tesa) — = eV? expf —e0 I), wo < y <0 
— co <p <oco (A.1) 

eg | 
Fly; p, ¢) = 1 — exp} —e Hin (4.2) 


where o = 1/8 ,u =1n @ and In is the natural logarithm. 
Letting L* denote the approximated likelihood function and 
Xp =(Y% — w/o, from Chapter II the approximated likelihood equations 


for wu and o are: 


dinL . dInL” | 


n—-S 
s(1 — On—s Br—sX(n—s)) ve (a; a iA = 


36 


dInL _ 6InL” 


Be we SE A + Xgl — 5X41) - (4.4) 


oo An—s + Pn—s(n=s) ) = ha (a; — B; iA(y)] = = 


It is not possible to compute the conditional s-bias of @ exactly nor 
the conditional s-bias of ft. The conditional s-bias of ¢ however, can be 


approximated by : 


EL Siok ] 
(A.5) 
rc ing 1 


Working out this approximation is tedious, but Balakrishnan and 
Varadan [ Ref. 6 p-149 ] provide a table of constants for the conditional 
s-bias of ¢. 

Let C denote the bias value of o from the table. Then 


El¢]=at+o0C =(14+Co. 


Thus o* = a/(1 + C) is nearly unbiased estimator of a. 


a 


APPENDIX B. DERIVATION OF EQUATION FOR DEGREES OF 
FREEDOM 


One basic expression used to determine the degrees of freedom, DF, 
for the chi-square distribution when obtaining confidence limits on the 
reliability of a complex system is DF = 2NFC, where NFC is total 
number of failed components. 

In chapter V, the analysis of the simulation results show that another 
formulae for computing degrees of freedom yielded a more accurate lower 
confidence limit procedure. The formula DF = 2NFC - NCOMP, for 
example, gave more accurate results than procedures that used DF = 
2NFC. Moreover, the accuracy associated with any fixed formula for DF 
degraded as the extent of truncation increased. That is, if for a particular 
component, NC items are tested until NF fail, the accuracy of each pro- 
cedure decreased as NF decreased. The decrease in accuracy became sig- 
nificant for values of NF < NC/2. Consequently, a formula for DF was 
developed that included NC/NF as one of the terms in the expression for 
DE: 

This formula for DF has the form DF = 2NFC + c(NC/NF) where 
NF is the smallest failure truncation value for each component type, NC 
is the number of items placed on test for that component with smallest 
NF and c is an unknown constant. An equation was established using this 
formula for each of the series system cases simulated using the value of 
DF that yielded the most accurate result and the appropriate values of 
NFC, NC and NF for that case. Each equation can be solved for c. Since 


there were thirty series system cases simulated (See Table 1A through 
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Table 3B), the resulting thirty equations yielded thirty values of c. An 


averaging process using least square methods was used to determine one 


value for c. The following set of equations show these computations for 


each confidence level value. 


A. ALPHA = 0.1 (CONFIDENCE 90 % ) 


Sy 
Zoe 
168 
104 

40 
eS 2 
p32 
168 
104 
40 
144 
224 
160 
96 

64 
144 
224 
160 
lez 

64 
160 
256 


60 + c¢ reduces ters = c 

240 + c reduces to -8 = c 

176 + I36c¢ reduces to:-8 = |.36c 
112 + 2.14c reduces to -8 = 2.14c 
45e+ S.0c reduces to -8 = 5.0c 
160 + c reduces to -8 
240 + c reduces to -8 
176 + 1.36c reduces to -8 l36c 
112 + 2.14c reduces to -8 = 2.14c 
48 + 5.0c reduces to -8 = 5.0c 
IGOes semcduces 10 -16 —.c 

240 + c reduces to -16 = c 

176 + 1.36c reduces to -16 = 1.36c 
112 + 2.14c¢ reduces to -16 = 2.14c 
48 + 5.0c reduces to 16 = 5.0c 
160 + c reduces to -16 
240 + c reduces to -16 = c 

176 + 1.36c reduces to -16 = 1.36c 
[12 ty2altc reducesito 0 = 2.14c 
48 + 5.0c reduces to 16 = 5.0c 
160 + c reduces to0 = c 

240 + c reducesto 16 =c 


Cc 


Cc 


C 


a7 


192 = 176 + 1.36c¢ reduces to 16 = 1.36c 
128 = 112 + 2.14c reduces to 16 = 2.14c 
64 = 48 + 5.0c reduces to 16 = 5.0c 

152 = 160 + c reduces to -8 
232 = 240 + c reduces to -3 
176 = 176 + 1.36¢ reduces to 0 = 1.36c¢ 
112 = 112 + 2.14c¢ reduces to 0 = 2.14c 
48 = 48 + 5.0c reduces to 0 = 5.0c 


The above equations form a linear system that is overdetermined. By 


C 


| 


C 


using the normal equation, A’4 x = A’b, we then find the least square 
solution to that overdetermined system which yields c = 1.1. Replacing 
the constants in the general formula, we have the approximating equation 
for obtaining the degrees of freedom ; 

DF = 2NFC + 1.1(NC/NF) 


B. ALPHA = 0.2 (CONFIDENCE 80 % ) 


152 = 160 + c reduces to -8 =c 
232 = 240 + c reduces to -8 = Cc 
168 = 176 + 1.36¢ reduces to -8 = 1:36¢ 
104 = 112 + 2.14c reduces to -8 = 2.14c 


40 = 48 + 5.0c reduces to -8 = 5.0c 


152 = 160 + c reduces to -8 =c 
232 = 240 + c reduces to -8 = ¢ 
168 = 176 + 1.36c reduces to -8 = 1.36c 
104 = 112 + 2.14c reduces to -8 = 2.14c 


40 = 48 + 5.0c reduces to -8 = 5.0c 
144 160 + c reduces to -16 = c 
224 = 240 + c reduces to-l6=c 


| 
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160 = 176 + 1.36c reduces to -16 = 1.36c 
96 = 112 + 2.14c reduces to -16 = 2.14c 
64 = 48 + 5.0c reduces to 16 = 5.0c 
144 = 160 + c reduces to -16 
224 = 240 + c reduces to -16 
160 = 176 + 1.36c reduces to -16 = 1.36c 
96 = 112 + 2.14¢ reduces to -16 = 2.14c 
64 = 48 + 5.0c reduces to 16 = 5.0c 
Ip2 — 160" c.réducés to -8 = c 

240 = 240 + creduces toO0 =c 

176 = 176 + 1.36c reduces to 0 | .36¢ 
2 = 112 + 2.14c reduces to 0 = 2.14c 
48 = 48 + 5.0c reduces to 0 = 5.0c 
l>z)— 160 4 ¢ reduces to -8 = c 

232 = 240 + c reduces to -8 = c 

168 = 176 * 1|.36c reduces to -8 = 1.36c 
2) = 112 3 2.l4e reduces to 0 = 2.)4c 
48 = 48 + 5.0c reduces to 0 = 5.0c 


Using the normal equation, A’4 x = A’b again, the least squares 


C 


C 


solution isc = 1.0. The approximate equation for obtaining the degrees 
OF freedom is DF = 2NFC + (NC/NF) 
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APPENDIX C. USERS’ GUIDE FOR RETP 


Reliability Estimation Test Plan (RETP). 





1. Brief Description 
RETP is a computer program written in FORTRAN that runs 
on the Amdahl mainframe at NPGS. It allows the user to simulate | 
exponential Weibull failure times of component items being tested 
to evaluate the accuracy of a confidence limit estimation procedure 
based on Type II data censoring (that is, testing 7, items of compo- 
nent 1 until f, of them fail). 


2. Program Input. (INI.DATA) 
The inputs of the program are specified to the program via an 
input file called INI.DAT. A sample input file 1s shown below. 
This sample input refers to a system with configuration as re- 
presented by Figure 4.4. 


SS SSS SS SS SS Se SSS SS 


This file contains the inputs required by the RETP model. 
Update only the numerical values between dotted lines as appropriate. 
Do not delete any of the comment lines. (INI. DAT) 


3 INT - Configuration of the system CONFIG 
1 = all subsystems in series 


2 = all subsystems in parallel 
3 = subsystems are connected 
in series-parallel 

16807.0 REAL = Initial random seed ISEED 
15 INT - Total # of subsystems in system NCOMP 
10 INT - # of exponential subsystems NEXP 
5) INT = ## of Weibull subsystems NWEI 
0 INT = ## of geometric subsystems NGEO 
0.10 REAL - Desired significance level ALPHA 
1000 INT = # of replications desired NREP 
3 INT = Test case number TCN 

1 = all exponential 

2 = all Weibull 

33= Eee + WEI 

4, = EXP + WEI + GEO 
15 INT = Number of cut sets NCS 
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c TEST PLAN : Testing NC(I) items of component i until NF(i) of them 
Cc roils 


c (Use REAL numbers ONLY !!!) 
rere. ee rere: 2 eee ee SS 8  e-- 
C Components Configuration 
c CSERIES( i) CSSGROUP( i) CSGROUP(i) NCGROUP( i) 
re er ee eee 8 ei ii we eee eee 
2 4 8 14 
Zz 4 8 14 
Z. 4 8 14 
Z 4 8 14 
Zz 4 8 14 
Zz 4 8 14 
2 4 8 14 
Z 4 8 14 
Z. 4 6 14 
2 4 6 14 
Zz 4 6 14 
2 4 6 14 
1 a 6 14 
1 1 6 14 
1 1 i 1 
i er Ree ee i | on eo ie ee - oe - - ~~~ 
c Comp Comp Comp Parameters ies tT Test Plan 
c Number Type Scale Shape Time/Cycle input 
ao TL ice.) PARM(1,1I) PARM(2,1) Ul@m § Ue(I) NC(I) NFCI) 
eeeint Int Real Real (hrs) Int Int Int 
es § Ser 35 s55eg55 49 8 ee a 
Jas 2.0 0.07500 hg TOn0 ZO; 0 15.0 
20 Ze 0.05500 Zao 8.0 2070 15.0 
3.0 220 0.17500 Je 8.0 20. 0 pels..0 
4.0 Ze) UrersU0 ono ga Z050— 15..0 
oy 18 20 UT IZs00 20 Sa0 ZOU) 15,0 
6.0 10 0.05500 ilegf0 0 20.0 1520 
i220 10 0.15000 1.0 TO ZnO l5e,0 
8.0 i are 0. 15500 130 2.0 ZO a) 155 0 
210 0 0.07500 10 20) ZO 1520 
10.0 ewe, 0.05000 30 230 Z0m0) 15,0 
M10, ie 0.17500 0 ZO Z0me 5" 0 
PZ 0 1550 0. 08500 1.0 lie, > 200m 1580 
digs). 10 0 0.17500 aU ro Z0—) 15.0 
14.0 150 0.05000 leno 1.0 2030515. 0 
>. 'U 0 0.02500 1330 bags 700m >. 0 
le ee ee Se. Se ee eee 
c Note : TY(I)=1 EXPONENTIAL P(surv) = exp(-PARM(1,I) * T) 
S TY(I)=2 WEIBULL P(surv) = exp(-(PARM(1,1)*T)**PARM(2,I)) 
c TY(I)=3 GEOMETRIC P(surv) = PARM(1,1) ** T 
i 8 
c SYSTEM CONFIGURATION : Identification of CUT SETS 
Cc Minimum groups of subsystems that have to fail 
Cc for the system to fail. 
ra a ale ae Se wc oo En wn nw SSS wo BE eee ewe ane 
ec Cut Set # in Set List of Components in Cutset 
cas Covel oly seOMr(y,Z2) ..... up to COMP(J,1) components 
eee = ose ee ee she 6 oe eee a ee 5 eee 
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‘®) 


oO 


oanananaoanananananaananaanaoanaaaan 


1 u i 0 0 0 0 0 0 0 0 
Z i Z 0 0 0 0 0 0 0 0 
S I 3 0 0 0 0 0 0 0 0 
4 1 4 0 0 0 0 0 0 0 0 
5 1 5 0 0 0 0 0 0 0 0 
6 1 6 0 0 0 0 0 0 0 0 
7 1 ii 0 0 0 0 0 0 0 0 
8 1 8 0 0 0 0 0 0 0 0 
9 1 9 0 0 0 0 0 0 0 0 
10 i 10 0 0 0 0 0 0 0 0 
liga ih 11 0 0 0 0 0 0 0 0 
2 1 WZ 0 0 0 0 0 0 0 0 
13 1 nS) 0 0 0 0 0 0 0 0 
14 1 14 0 0 0 0 0 0 0 0 
es 1 tS 0 0 0 0 0 0 0 0 


3. Program Flow and _ Logic. (NAME.DEF,PARM.DEF and 
RETP.FOR) 

Input parameters are first read in by the program by calling the 
INPUT subroutine. The program then evokes the SIM subroutine 
which generates the random failure times and computes the key sta- 
tistics required in the procedure. The next subroutine EVAL deter- 
mines the measures of accuracy for each case. REPORT is the 
subroutine which generates the output file for the run OUT.DAT. 

The variables in the program RETP.FOR are described in the file 
NAME.DEF as listed below. 


This file contains the declaration for input and output variables 
used in the RETP model. (NAME. DEF) 


Input Variables. 

ISEED = initial random seed selected. (Real) 

SEED = current random seed (Real) 
RS = true overall system reliability. 
ALPHA = desired significance level (Real) 
NREP = number of replication desired for simulation (Integer) 
TPN = test plan number 
TCN = test case number (1, 2, 3 or 4) 
NCOMP = total # of components in the system (Integer) 
NEXP = number of components with EXP failure times. (Integer) 
NWEI = number of components with WEI failure times. (Integer) 
NGEO = number of components with GEO failure times. (Integer) 
Distribution: EXPonential WEibull GEOmetric 
TY(Ci) = type: 1 2 3 

PARMO@ ie). ©: Scale(1/hr) Scale(1/hr) Prob 
PARM(2,i) - Shape = 
UTC i) = utilization time (hrs) for component i (EXP and WEI) 
UC(i) = utilization cycles for component i (GEO only) 
NC( i) = number of test samples (sample size) for component i. 
NF(i) = desired number of failures in test for component i. 


oononnoannanananannnnn 


OQ 


oaoonoanoanannnnnnnnnnnannn 


NCS 
COMPC I,K) 


CONFIG 
CSERIES(i) 
CSGROUP( i) 


CSSGROUP(i) 


NCGROUP( i) 


number of cut sets for the system. 

kth parameter of cutset j ( first being the number of 
components belonging to the cut-set) 

configuration number of components arrangement( integer). 
number of subsystems in series in a subsubgroup(integer) 
number of subsystems i in a subgroup (integer) 

number of subsystems i in a subsubgroup (integer) 

number of subsystems i in a group (integer) 


Assumed Variables. 


MAXCOMP = maximum number of components allowed in the system 

MAXREP = maximum number of replication permitted. 

MAXCUT = maximum number of cut-sets. 

Program and Output Variables. 

TT(i) = total accumulated failure time (in hour) for component i 
(EXP and WEI only) 

TC(i) = total accumulated cycles to failure (incl. failure cycle) 
for component i (GEO only) 

EBETA(i) = estimate for shape parameter of component i (if Weibull) 

REL1(j) = actual reliability for cut-set j. 

REL2( j) = computed reliability for cut-set j for current replication 

ELM(i) = estimated component failure rate (li/hrs) for component i 

ELMAX(m) = max. estimated component failure rate for rep. m (1/hrs) 

ERC i) = ratio of estimated failure rate to ELMAX. 

NFC(m) = total number of failed test components. 

LMU(m) = upper confidence limit for failure rate (1/hrs). 

RSL(m) = lower confidence limit estimated for system reliability 
for the mth replication. 

ORSL(m) = ascending order of RSL(M). 

RSLOW = (1-ALPHA)x100 percentile of set of RSL(M). 

LEVEL = achieved confidence interval i.e proportion of RSL(M) 
that are lesser than RS (conservative estimate) 

-- END OF NAME. DEF creer ccr reser ce reer cence ccc en ccc ener ens ernnsensscan- 


Together with the main program in RETP.FOR are the other sub- 


routines needed in the simulation. The declaration of variables is done in 


the file PARM.DEF. Relevant descriptions are included as comment lines 


in the source code to help explain the program segments. A listing of 
a. DEF and RETP.FOR is given below. 


Cc 


This file contains the declaration for input and output variables 
c¢ used in the RETP model. ( PARM. DEF) 


INTEGER MAXCOMP, MAXREP 
PARAMETER( MAXCOMP = 100, MAXREP = 1000, MAXCUT = 20 ) 
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REAL*8 ISEED, SEED 

INTEGER NREP, TCN, NCOMP, NEXP, NWEI, NGEO, NCS, 
NC(MAXGOMP) ,NF(MAXCOMP), TY(MAXCOMP), NFC(MAXREP), 
UC(MAXCOMP), TC(MAXCOMP), COMP(MAXCUT,MAXCOMP) , 
CONFIG, CSERIES(MAXCOMP), CSSGROUP(MAXCOMP) , 
CSGROUP( MAXCOMP) , NCGROUP( MAXCOMP) , DFR 

REAL*8 RS, ALPHA, UT(MAXCOMP), TT(MAXCOMP), 
PARM(2,MAXCOMP), ELM(MAXCOMP), ER(MAXCOMP), 
LMU(MAXREP), RSL(MAXREP), ORSL(MAXREP), 
ELMAX(MAXREP), RSLOW, LEVEL, EBETA(MAXCOMP) , 
REL1(MAXCUT), REL2(MAXCUT) 


t++++ 


t++++ 


{COMMON/ BLOCK1/1SEED, SEED, NREP, TCN, NCOMP, NC, NF, NEXP, NWEI, 
NGEO, NCS, TY, NFC, UC, TC, COMP, CONFIG, CSERIES, 
CSSGROUP, CSGROUP, NCGROUP, DFR 

COMMON BLOCK2/RS ALPHA, UT, TT, PARM, ELM, ER, LMU, 
RSL, ORSL, ELMAX, RSLOW, LEVEL, EBETA, REL1, REL2 


+ + 


This file contains the main program and the subroutines 
for the Reliability Estimation Test Plan (RETP) model. 


IBM Mainframe version. 
Test Plan 1: Testing NC(I) items for component i 


PROGRAM RETP 
Include the declaration files. 


INCLUDE ‘NAME DEF' 
INCLUDE 'PARM DEF' 


Read the input data. 
CALL INPUT 
Activate simulation 
CALL SIM 
Process and evaluate output data. 
CALL EVAL 
Generate simulation report. 
CALL REPORT 


STOP 
END 
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c 2. Input Initialisation Subroutine (INPUT). 


SUBROUTINE INPUT 
This subroutine reads in the inputs for the RETP model. 
Include the declaration file. 

INCLUDE 'PARM DEF’ 


INTEGER I, J, K, DUM2(11), DUM3(4) 
REAL*8 DUM1(7) 


Read data from file 'INI.DAT' designated as logic unit 1. 
OPEN (UNIT=1,FILE='/INI data Al’) 


READ(1,*) CONFIG 
READ(1,*) ISEED 
READ(1,*) NCOMP 
READ(1,*) NEXP 
READ(1,*) NWEI 
READ(1,*) NGEO 
READ(1,*) ALPHA 
READ(1,*) NREP 
READ(1,*) TCN 
READ(1,*) NCS 


DO 40 I = 1, NCOMP 
READ(1,*) DUM3 


CSERIES(I) = DUM3(1) 

CSSGROUP(I) = DUM3(2) 

CSGROUP(I) = DUM3(3) 

NCGROUP(I) = DUM3(4) 
40 CONTINUE 


DO 50 K = 1, NCOMP 
READ(1,*) DUM1 


I = NINT(DUM1(1)) 
TY(I) = NINTCDUM1(2)) 


IF (TY(I).EQ.1) THEN 
PARM(1,1) = DUM1(3) 
PARM(2,1) = DUM1(4) 


UT(I) = DUM1(5) 

NC(I) = NINTCDUM1(6)) 
NF(I) = NINTCDUM1(7)) 
EBETACI) = 0 


ELSEIF (TY(1).EQ.2) THEN 
PARM(1,I) = DUM1(3) 
PARM(2,1) = DUM1(4) 


UT(I) = DUM1(5) 
NC(I) = NINT(DUM1(6)) 
NF(I) = NINT(DUM1(7)) 


ELSEIF (TY(I).EQ.3) THEN 
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PARM(1,1) = DUM1(3) 

PARM(2,1) = DUM1(4) 

UC(I) = NINT(DUM1(5)) 

UT(I) = DUM1(5) 

NC(I) = NINT(DUM1(6)) 

NF(I) = NINT(DUM1(7)) 
ENDIF 


50 CONTINUE 


DO 80 I = 1, NCS 
READ(1,*) DUM2 
J = DUM2(1) 
COMP(J,1) = DUM2(2) 
DO 70 K=1,COMP(J,1) 
COMP(J,K+1) = DUM2(K+2) 
70 CONTINUE 
80 CONTINUE 
CLOSE( UNIT=1) 
RETURN 
END 


SUBROUTINE SIM 


This subroutine simulates NREP possible outcomes of the test plan 
desired in order to obtain the raw estimates of LMU(M) and RSL(M) 
for each of the replication. 


Include the declaration file 
and declare local variables. 


INCLUDE 'PARM DEF' 


INTEGER I, J, K, M, ISUM, KEY, L, C, CC, ccc 
REAL*8 UNI 

REAL*8 SUM, PROD, RP, RPP, RSP, RSS, 
i FT(MAXCOMP), OFT(MAXCOMP) , YY(MAXCOMP) 


LOGICAL TEMP 
SEED = ISEED 
Compute overall true system reliability RS. 


C= 0 

CC = 0 
CCC = 0 
RP = 0.0 
RPP = 1.0 
RSP = 0.0 
RSS 
RS 
J=1 

TEMP = . FALSE. 

DO WHELEC. NOT ALEHP ) 


hes, 
Ls, 
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DO 20 L= J, J + (CSERIES(J) - 1) 
PROD = 1.0 
DO 10 I = 1, COMP(L,1) 
K = COMP(L, I+1) 
PROD = PROD*(1- SURV(TY(K) ,PARM(1,K),PARM(2,K),UT(K))) 
10 CONTINUE 
RELI(L) = 1.0 - PROD 
RSS = RSS * RELI(L) 


20 CONTINUE 


L=L-1 
RP = 1.0 - (1.0 -RSS)*(1.0 -RP) 
RSS =1.0 
IF (L.GE.(C + CSSGROUP(L))) THEN 
RPP = RPP * RP 
RP = 0.0 
C = C + CSSGROUP(L) 
ENDIF 
IF (L.GE.(CC + CSGROUP(L))) THEN 
RSP = 1.0 - (1.0 - RPP)*(1.0 - RSP) 
RPP = 1.0 


CC = CC + CSGROUP(L) 
ENDIF 


IF (L. GE. (CCC + NCGROUP(L))) THEN 
RS = RS * RSP 


RSP = 0.0 
CCC = CCC + NCGROUP(L) 
ENDIF 


J = Je+eCSERIESCL) 
IF (L. GE. NCOMP) TEMP = . TRUE. 
ENDDO 


Start of Simulation 
(Initialize replication counter M) 


M=1 
DO WHILE (M. LE. NREP) 


Test Plan : Sample and determine unknown TT(I) 
ceccteeee with known NC(I) until NFCI) fails. 


Generate NC(I) failure times, put them in ascending order 
with the smallest failure time on the top of the list. 


DO 70 I = 1, NCOMP 


DO 40 K = 1, NC(I) 
CALL LRNDPC(SEED,UNI, 1) 
IF (TY(I).EQ. 1) THEN 
FT(K)= -LOG(UNI)/PARM(1,1) 
ELSEIF (TY(1I).EQ.2) THEN 
FT(K)=(1. 0/PARM( 1,1) )*( -LOGCUNI) )**(1. 0/PARM(2,I)) 
ELSEIF (TY(1).EQ.3) THEN 
FT(K)= 1.0 
DO WHILE (UNI. LT. PARM(1,1)) 
FT(K) = FT(K) + 1.0 
CALL LRNDPC(SEED,UNI,1) 
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ENDDO 
ENDIF 
CONTINUE 


Buble Sort the failure times in ascending order. 


CALL BUBBLE(NC(I),FT,OFT) 


Take logarithm of the ordered failure times. 


DO 45 K=1, NCCI) 


YYCK) = LOGCOFT(K)) 


45 CONTINUE 


Compute the total time accumulated in the test and the estimate 
for the failure rate of the component as in the procedure. 
eens (1) NE. 2) THEN 
SUM = 0.0 
DO 50 K = 1, NF(I) 
SUM = SUM + OFT(K) 
a0 CONTINUE 
TTC I) = FLOATCNC(I)-NFCI))*OFTCNF(I)) + SUM 
TP Cty ( 1). EO.1)- aHEN 
ELM(I) = FLOAT(NF(I) - 1)/TTCI) 
ELSEIF (TY( Ty} HOs3)° THEN 
ELM(I) = FLOATCNF(I) - 1)/TTCI) 
ENDIF 
EBSEIP (TYC1). EQ.2) THEN 
CALL APROXMLE(CYY,NC( I) ,NFCI) ,EBETACI)) 
EBETAC(I) = BNCNC(I))*EBETACT) 
SUM = 0.0 
DO 60 K = 1, NF(CTI) 
SUM = SUM + OFT(K)**EBETA(T) 
60 CONTINUE 
TT(I) = FLOATCNC(I)-NFCI))*OFTCNF(I))**EBETACI) + SUM 
ELM(I) = FLOATC(NFC(1I))/TTCI) 
ENDIF 
70 CONTINUE 
Determine the total number of failed test items. 
ISUM = 0 
DO 80 I = 1, NCOMP 
ISUM = ISUM + NF(T) 
80 CONTINUE 


NFC(M) = ISUM 


Determine the maximum failure rate estimate 


70 
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Cc 


and identify the component. 


90 


ELMAX(M) = 0.0 
KEY = 0 
DO 90 I = 1, NCOMP 
IF (ELM(1). GT. ELMAX(M)) THEN 
ELMAX(M) = ELM(I) 
KEY = I 
ENDIF 
CONTINUE 


Compute the ratios of the failure rate estimate to their maximum 


100 


DO 100 I = 1, NCOMP 
ER(I) = ELM(I)/ELMAX(M) 
CONTINUE 


c Determine LMU(M) 


Cc 


Cc 


110 


SUM = 0.0 
DO 110 I = 1, NCOMP 

SUM = SUM + (ER(I)*TT(I)) 
CONTINUE 


DFR = 2*NFC(M)-2*NCOMP 
LMU(M) = CHISQD(1-ALPHA, DFR) /(2*SUM) 


c Compute estimate of overall reliability RSL(M) for the system. 


c 


b15 


Ale) 


c= 0 
CCl=— 20 
cole = 0 
RP = 0.0 
RPP = 1.0 
RSP = 0.0 
RSS= 1.0 
RSL(M) = 1.0 
i= fi 
TEMP = . FALSE. 
DO WHILE(. NOT. TEMP) 
DO 120 L = J, J + (CSERIES(J) - 1) 
PROD = 1.0 
Doelis. 1 = 1, COMP( 1) 
K = COMP(L,I+1) 
IF (TY(K).EQ.1) THEN 
PROD = PROD*(1-SURV(TY(K) , LMU(M)*ER(K) , EBETA(K) , UT(K))) 
ELSEIF (TY(K).EQ. 2) THEN 
PROD = PROD*(1-SURV(TY(K) , (LMU(M)*ER(K))**(1. /EBETA(K)), 


re EBETA(K) ,UT(K))) 


FlSruemcarck). EQ.3) THEN 
PROD=PROD*( 1-SURV(TY(K) , 1. DO-LMU(M)*ER(K) ,0. DO,UT(K))) 
ENDIF 
CONTINUE 
REL2(L) = 1.0 - PROD 
RSS = RSS*REL2(L) 
CONTINUE 


mM 
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L = 1L =a 
RP = 1.0 - (1.0 - RSS)*(1.0 - RP) 
RSS =1.0 
IF (L.GE.(C + CSSGROUP(L))) THEN 
RPP = RPP * RP 
RP = 0.0 
C = C + CSSGROUP(L) 
ENDIF 
IF (L.GE. (CC + CSGROUP(L))) THEN 
RSP = 1.0 - (1.0 - RPP) * (1.0 - RSP ) 


RPP = 1.0 
CC = CC + CSGROUP(L) 
ENDIF 


IF (L.GE. (CCC + NCGROUP(L))) THEN 
RSL(M) = RSL(M) * RSP 


RSP = 0.0 
ccc = CCC + NCGROUP(L) 
ENDIF 


J = J + CSERTIESCL) 
IF (L.GE.NCOMP) TEMP = . TRUE. 
ENDDO 
Increment replication counter. 
M=M+t+1 
ENDDO 


RETURN 
END 


SUBROUTINE EVAL 


This subroutine calls BUBBLE to sort the array RSL(NREP) in 
ascending order to get an ordered array ORSL(NREP). It also 
determine the estimate for RSLOW at the specified significance 
level ALPHA and the value of Level in which ORSL(LEVEL) is closest 
to the true reliability RS. 


Include the declaration files 
and declare the local variables. 


INCLUDE 'PARM DEF' 


INTEGER INDEX, M 
REAL*8 DIFF 


Order the array RSL(NREP) in ascending order. 


DO 10 M = 1, NREP 
ORSL(M) = RSL(M) 


10 CONTINUE 


12 


'?) 


Bubble Sort. Sink the larger of the pair. 


'?) 


CALL BUBBLE(NREP,RSL,ORSL) 


‘a 
c Determine the (1-ALPHA) % lower confidence bound for the system 
c reliability. 
‘a 
RSLOW = ORSLC(NINTCNREP*( 1-ALPHA) )) 
C 
c Finding the % confidence level for the true reliability RS. 
c (i.e the proportion of RSL(M) lesser than RS) 
c 
DIFF = 1.0 
INDEX = 0 


WOEzZOO M =.1.. NREP 
IF (ABSCORSL(M) - RS). LT.DIFF) THEN 
DIFF = ABSCORSL(M) - RS) 
INDEX = M 
ENDIF 
200 CONTINUE 


LEVEL = FLOAT( INDEX) /NREP 
a 


2) 


Record evaluated parameters in RAW1.DAT (unit 2). 


OPEN( UNIT=2) 
WRITE( 2,300) 
300 FORMAT(1x, ' M LMU(M) ELMAX(M) RSL(M)', 
a ORSL(M) NFC(M)') 
DO 500 M = 1, NREP 
WRITE(2,400) M,LMU(M) ,ELMAX(M) ,RSL(M) , ORSL(M) ,NFC(M) 
400 FORMAT(1x,16,2F12.7,2F12.7,16) 
500 CONTINUE 
CLOSE(UNIT=2) 


RETURN 
END 
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SUBROUTINE REPORT 


oO 


This subroutine record the simulation result into the 'OUT DATA’ 
file indicated by logic 3. 


Include the declaration files 
and declare local variables. 


eooannnno 


INCLUDE ' PARM DEF’ 
INTEGER I, J, K, DUM(10) 


Write to output file 'OUT DATA’ designated as logic unit 3. 


oO oO 


OPEN( UNIT=3) 
WRITE(3,10) 
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200 


300 


400 


500 


10 
20 
Zo 
26 
Za 
28 
2g 
30 
40 
50 
60 
70 
80 
90 


100 
110 
1Z0 


130 
140 


+ 


WRITE(3,20) NREP 

WRITE(3,25) NCOMP,CONFIG 

IF (CONFIG. EQ. 1) THEN 
WRITE( 3,26) 

ELSEIF (CONFIG. EQ. 2) THEN 
WRITE(3, 27) 

ELSE 
WRITE(3, 28) 

ENDIF 

WRITE(3,29) DFR 

WRITE(3, 30) 

WRITE(3,40) 

WRITE(3,50) ISEED,NCOMP,ALPHA,NCS,TCN 


WRITE(3,60) 
DO 200 I = 1, NCOMP 

WRITE(3,70) I,TY(I),PARM(1,1),PARM(2,1),UT(1I),NC(1I) ,NF(I) 
CONTINUE 
WRITE(3,80) 
WRITE( 3,90) 
DO 300 I = 1, NCOMP 

WRITE(3,100) I,NF(I),TT(1I),ELM(I),ER(I),EBETA(I) 
CONTINUE 
WRITE(3,110) 
WRITE(3,120) 
DO 500 J = 1, NCS 

DO 400 K = 1, 10 

DUM(K) = COMP(J,K) 

CONTINUE 
WRITE(3,130) J,DUM,REL1(J) ,REL2(J) 
CONTINUE 
WRITE(3,140) 
WRITE(3,150) RS,ELMAX(NREP) , LMU(NREP) , RSLOW, LEVEL 


FORMAT(1x,'OUT1 DATA : Output File of the RETP1 simulation’ ) 
FORMAT( 1x, ' after ',15, repWiation ,/) 
FORMAT(1x,'COMMENTS : ',12,' COMPONENTS IN CONFIGURATION: ',1I2 ) 
FORMAT( 1x, ' (SERIES SYSTEM) /.) 

FORMAT(1x,' (PARALLEL SYSTEM) ee) 

FORMAT( 1x, ' (SERIES-PARALLEL SYSTEM) oe) 
FORMAT(1x,' DF = 2 * (NFC - NCOMP) =',I4 ) 
FORMAT( 1x, ‘Input parameter: ',/) 

FORMAT(1x,' ISEED NCOMP ALPHA NGS TCM’ WW) 
FORMAT(1x,F10.1,18,F8.4,216,/) 

FORMAT(1x,’ I. TY(1I) PARM1(I) PARM2(I) UT(I) NC(L) NE@L) ae 
FORMAT( 1x,12,16,1x,2F9.5,F8.2,2x,216) 

FORMAT(1x,/, Output Parameters for the Last Replication: ',/) 
FORMAT(1x,' I NF(I) Tien) ELM(I) ERG 


in '  EBETA(I)',/) 


FORMAT(1x,12,16,4x,E14.7,2x,F9. 7,2x,F9. 7,2x,F9. 7) 
FORMAT(1x,/,'Cut-Set Data: ',/) 
FORMAT(1x,' J NUM Component List - 


+ ' — REL1(J) RELZCM) = 6=6— oe 


FORMAT(1x,1I2,I5,913,2F12.9) 
FORMAT(1x,/,' RS ELMAX(M) LMU(M) ', 
RSLOW LEVEL ',/) 
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P50 “POURMAT Cie oF i2.7,/) 
CLOSECUNIT = 3) 


RETURN 

END 
This portion of the file contains functions and subroutines 
used in the RETP model. 


SUBROUTINE LRNDPC (DSEED,U,N) 


INTEGER N, I 
REAL*8  U(N) 
REAL*8 D31M1, DSEED, D31 


DATA D31M1/2147483647. DO/ 
DATA D31 /2147483648. DO/ 
DO 5. I = 1,!N 
DSEED = DMOD( 16807. DO*DSEED ,D31M1) 
5 U(I) = DSEED / D31 
RETURN 
END 


FUNCTION SURV( TYPE, PAR1,PAR2,UTIL) 


This function returns the survival probability of the component of 
different types (TYPE) with scale (PAR1) and shape (PAR2) parameters 
given the specified utilization times or cycles (UTIL). 


eyercin TYPE, N 
REAL*8  PAR1, PAR2, UTIL 


PCr Oe)) MEN 
SURV = EXP(-( PAR1*UTIL) ) 
Piehif (CTYPE. EQ. 2) THEN 
SURV = EXP(-((PAR1*UTIL)**PAR2) ) 
ELSE 
N = NINTCUTIL) 
SURV = PAR1**N 
ENDIF 
END 


SUBROUTINE BUBBLE(N,LIST,OLIST) 
This subroutine performs a bubble sort in increasing order. 


(i.e sink the greater numeral) for the first N terms in an array 
LIST and returns the result in OLIST. 
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LOGICAL DONE 

INTEGER N, K, PAIR 

REAL TEMP 

REAL*8 LIST(*), OLIST(*) 


Sink the larger of the pair. 


pO 50 K=1, N 
OLIST(K) = LIST(K) 
50 CONTINUE 
PAIR=N - 1 
DONE = . FALSE. 
DO WHILE (.NOT. DONE) 
DONE = . TRUE. 
DO 100 K = 1, PAIR 
IF (OLIST(K). GT. OLIST(K+1)) THEN 
TEMP = OLIST(K) 
OLIST(K) = OLIST(K+1) 
OLIST(K+1) = TEMP 
DONE = . FALSE. 
ENDIF 
100 CONTINUE 
PAIR = PAIR - 1 
ENDDO 
RETURN 
END 


FUNCTION BN(T) 


This function returns the value of the unbiased factor for the biased 
approximate MLE of the shape parameter of a Weibull distribution 
with a sample size of N. 


INTEGER I 


IF (I. LE.5) THEN 
BN = (I*0.699)/(5.0) 

ELSEIF (I.EQ.6) THEN 
BN = 0.752 

ELSEIF (I.EQ.7) THEN 
BN = 0.786 

ELSEIF (I.EQ.8) THEN 
BN = 0.82 

ELSEIF (I.EQ.9) THEN 
BN = 0.8395 


ELSEIF (1. EQ. 10) THEN 
BN = 0.859 

ELSEIF Cl. EQ7i))) THEN 
BN = 0.871 

BESEIF (1.EQ.12) THEN 
BN = 0. 883 
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ELSEIF (I.EQ.13) THEN 

BN = 0. 892 
ELSEIF (I.EQ.14) THEN 

BN = 0.901 
ELSEIF (I.EQ.15) THEN 

BN = 0.9075 
ELSEIF (I.EQ.16) THEN 

BN = 0.914 
ELSEIF (I.EQ.17) THEN 

BN = 0.9185 
ELSEIF (I.EQ.18) THEN 

BN = 0.923 
ELSEIF (I.EQ.19) THEN 

BN = 0.927 
ELSEIF (I.EQ.20) THEN 

BN = 0.931 
ELSEIF (I. LE. 25) THEN 

BN = 0.931 + (I-20)*0.014/5.0 
ELSEIF (I. LE. 30) THEN 

BN = 0.945 + (I-25)*0. 01/5. 0 
ELSEIF (I. LE.40) THEN 

BN = 0.955 + (I-30)*0. 011/10. 0 
ELSEIF (I. LE.60) THEN 

BN = 0.966 + (I-40)*0. 012/20. 0 
ELSEIF (I. LE. 80) THEN 

BN = 0.978 + (I-60)*0. 006/20. 0 
ELSEIF (I. LE. 100) THEN 

BN = 0.984 + (I-80)*0. 003/20. 0 
ELSEIF (I. LE. 120) THEN 

BN = 0.987 + (I-100)*0. 003/20. 0 
ELSE 

BN = 1.0 
ENDIF 
RETURN 
END 


FUNCTION BIAS(N,S) 


This function returns the value of the unbiased factor for the biased 
approximate MLE of the scale parameter of a extreme value distribution 
with a sample size of N. 


INTEGER N, S$ 


IF (N.EQ.10) THEN 

IF (S.EQ.0) THEN 
BIAS = 0.9339 

ELSEIF (S.EQ.1) THEN 
BIAS = 0.9275 

ELSEIF (S.EQ.2) THEN 
BIAS = 0.9152 

ELSEIF (S.EQ.3) THEN 
BIAS = 0.9001 


a 


ELSEIF (S.EQ.4) THEN 
BIAS = 0. 8908 
ELSEIF (S.EQ. 8) THEN 
BIAS = 0. 8453 
ELSE 
BIAS = 0.7998 
ENDIF 
ELSEIF (N.EQ.15) THEN 
IF (S.EQ.0) THEN 
BIAS = 0.95025 
ELSEIF (S.EQ.1) THEN 
BIAS = 0.94715 
ELSEIF (S.EQ. 2) THEN 
BIAS = 0.94015 
ELSEIF (S.EQ.3) THEN 
BIAS = 0.9309 
ELSEIF (S.EQ.4) THEN 
BIAS = 0.92495 
ELSEIF (S.EQ. 8) THEN 
BIAS = 0. 89855 
ELSE 
BIAS = 0. 8718 
ENDIF 
ELSEIF (N.EQ. 20) THEN 
IF (S.EQ.0) THEN 
BIAS = 0. 9666 
ELSEIF (S.EQ.1) THEN 
BIAS = 0. 9668 
ELSEIF (S.EQ.2) THEN 
BIAS = 0. 9651 
ELSEIF (S.EQ.3) THEN 
BIAS = 0.9617 
ELSEIF (S.EQ.4) THEN 
BIAS = 0. 9591 
ELSEIF (S.EQ. 8) THEN 
BIAS = 0.9518 
ELSE 
BIAS = 0. 9438 
ENDIF 
ENDIF 
RETURN 
END 


SUBROUTINE APROXMLE(CY,NN,A,BETAHAT) 


INTEGER I, NN, R, S, A 
REAL*8 Y(*), P(20), Q(20), ALPA(20), BETA(20), 
+ BETAHAT , GAMMA , DELTA,B,C,D,E,M,BB,CC,DD,EE,MM, SIGMAHAT 


~~ wa 


oil 


NN - K = A 
DO 5 I=1, A 
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P(I) = REAL(I)/(NN + 1) 
Q(T) iem= PCT) 
ALPACI)=1.0 + LOG(Q(1I))*(1.0 - LOG(-LOG(Q(I)))) 
BETACI)= -LOG(Q(T)) 
5 CONTINUE 


GAMMA = -(Q(R+1)/P(R+1) )*LOG(Q(R+1) )*(1. 0-LOG(-LOG(Q(R+1)))) + 
ar (QCR+1) /(C PC R+1)**2) )*( (LOG QC R+1) ) )**2)*LOG( -LOG(Q(Rt1) )) 
DELTA = (Q(R+1)/P(R+1))*LOG(Q(R+1) )*(1. 0+(LOG(Q(R+1))/P(R+1))) 
BB = 0.0 
CC = 0.0 
DD = 0.0 
EE = 0.0 
MM = 0.0 
DO 10 I=R+1, NN-S 
CC = CC + ALPA(I) 
BB = BB + BETA(I)*Y(I) 
MM = MM + BETA(I) 


10 CONTINUE 


M = R*DELTA + S*BETA(NN-S) + MM 
B = (R*DELTA*Y(R+1) + S*BETA(NN-S)*Y(NN-S) + BB)/M 
C = (R*GAMMA - S*(1.0 - ALPA(NN-S)) + CC)/M 


HORI SS tai | SeNN=S 
DD = DD + ALPACI)*(YCI) - B) 
ieee BEEACL) CC YC Ll) = B)**2) 
15 CONTINUE 


D 
E 


R*GAMMA*(Y(R+1)-B) - S*(1.0-ALPA(NN-S))*(Y(NN-S)-B) + DD 
R*DELTA*((Y(R+1)-B)**2) + S*BETA(NN-S)*((Y(NN-S)-B)**2) + EE 


SIGMAHAT =(-D + SQRT(D**2 + 4. O*XE*FLOAT(A)))/(2. O*XFLOAT(A) ) 
SIGMAHAT =SIGMAHAT/BIAS(NN,S) 
BETAHAT = 1. 0/SIGMAHAT 


RETURN 
END 


FUNCTION CHISQD(P,N) 


Modified version of Algorithm 451 from Communication of the ACM 
August 1977 Vol. 16 No. 8. 


This function evaluates the quantile at the probability level P 
(left tail area) for the Chi-Square Distribution with 
N degrees of freedom. 


REAL*8 P 
REAL X 
INTEGER IF 
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20 


30 


DIMENSION C(21), A(19) 


DATA C/ 1.565326E-3, 
1. 060438E-3, 
-6. 950356E-3, 
-1. 323293E-2, 
2.277679E-2, 
-8. 986007E-3, 
-1.513904E-2, 
2.530010E-3, 
-1.450117E-3, 
5. 169654E-3, 
-1.153761E-2, 
1. 128186E-2, 
2.607083E-2, 
-0. 2237368, 
9. 780499E-5, 
-8.426812E-4, 
3. 125580E-3, 
-8. 553069E-3, 
1. 348028E-4, 
0.4713941, 
1. 0000886 / 


t++tttteetteeteetteteteeet 


DATA A/ 1.264616E-2, 
-1.425296E-2, 
1. 400483E-2, 
-5. 886090E-3, 
-1.091214E-2, 
-2. 304527E-2, 
3.135411E-3, 
-2.728484E-4, 
-9.699681E-3, 
1. 316872E-2, 
2.618914E-2, 
=0. 2222220) 
5.406674E-5, 
3.483789E-5, 
-7, 274761E-4, 
3. 292181E-3, 
-8.729713E-3, 
0.4714045, 


Lawl 


PE ONeill 0 20 30) 

CALL XFROMP(. 5*(1. -P),X,IF) 
CHISQD = X 

CHISQD = CHISQD*CHISQD 
RETURN 


t++t+eteeetteteteet+ett 


CHISQD 
RETURN 


-2.*LOG(1. -P) 


F = N 
Fl = 1./F 
CALL XFROMP(P,X,IF) 
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T= xX 
F2 = SQRT(F1)*T 
IF (N. GE. (2+INT(4. *ABS(T)))) GO TO 40 


CHISQD = (((((((C(1)*F2+C(2))*F2+C(3))*F2+C(4) )*F2 
+C(5) )*F2+C(6) )*F24+C(7))*F1+((((((C(8)+C(9)*F2)*F2 
+C(10))*F2+C(11))*F2+C( 12) )*F24+C(13))*F2+C(14)))*F1+ 
(((C(C(15)*F24+C( 16) )*F2+C(17))*F2+C(18))*F2 
+C(19))*F2+C( 20) )*F2+C(21) 


t++tt 


GO TO 50 


40 CHISQD = (((A(1)+A(2)*F2)*F1+(((A(3)+A(4)*F2)*F2 
+A(5))*F2+A(6)))*F1+(((((AC7 )+A(8)*F2)*F2+A(9) )*F2 
+A(10))*F2+A( 11) )*F2+A(12)))*F1+(((((A(13)*F2 
+A(14))*F2+A( 15) )*F2+A(16))*F2+A(17))*F2*F2 
+A(18))*F2+A(19) 


t+++t 


50 CHISQD = CHISQD*CHISQD*CHISQD*F 


RETURN 
END 


SUBROUTINE XFROMP(P,X, IFAULT) 
mesocrtimeAs 24 °J.k.STAT. SOC. C (1969) Vol.18. No.3. 
This subroutine compute the standard normal deviate X for the 


specified left tail area P. 


REAL*8 P 
DIMENSION A(5) 
DIMENSION CONNOR (17), HSTNGS(6) 


+- 


DATA CONNOR/ 8. 0327350124E-17, 
1. 4483264644E-15, 
.4668270103E-14, 
.9554295164E-13, 
.9477940136E-12, 
.3507027951E-11, 
, 0892221037E-9, 
. 3122532964E-8, 
.4503852223E-7, 
.4589169001E-6, 
.3227513228E-5, 
. 0683760684E-4, 
.5757575758E-4, 
, 6296296296E-3, 
. 3809523810E-2, 
a 
. 33333333333 / 


Het tteteeteteteteeeest 
COON FN FR FR Re Fe © UI WW PO 
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DATA RTHFPI / 1. 
DATA RRT2PI / 0. 
DATA TERMIN / 1. 


DATA HSTNGS / 


t+++4++4+ 


IFAULT = 1 


OOrFOO fS 


2533141373 / 
3989422804 / 
OE-11 / 


5155105 


. 802853, 
. 010328, 
. 432788, 
. 189269, 
. 001308 / 


IF (CP. BeOS) SOR SCE. Gh. 1.0 jer eo 


IFAULT = 0 


Get first approximation Xo to deviate by Hasting's formula. 


B =P 
IF (BeCl.. 0.5) ae 


F 
E 


- LOG(B) 
SQRT(F+F) 


eet. = 2B 


XO = -E + (CHSTNGS(3)*E+HSTNGS( 2) )*E+HSTNGS(1))/ 
+ (C CHSTNGS(6)*E+HSTNGS(5 ) )*E+HSTNGS( 4) )*E+1. 0) 


EF <¢ 
XO = 0.0 
PO = 
Xl = 
GO TO 7 


XO. LT. 0. 0)"GGawoO 1 


Find the area PO corresponding to XO 


Ll Y¥ = XO**2 
IF (XO) 
Y = -0. 5*Y 


GOr1TO 3 


(1) series approximation 


PO = CONNOR(1) 
DO 2L=2, 17 

2 PO = PO*Y + CONNOR(L) 
PO = (PO*Y+1. 0)*XO 
X1 = -(PO+RTHFPI )*EXP(-Y) 
PO = PO*RRT2PI + 0.5 
Commo 7 


(2) continued fraction approximation 


3 wmeH 1.07% 


A(2) = 1.0 
AC(3) = 1.0 
A(4) = Z+ 1.0 
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A(5) = 1.0 


W = 2.0 
DO6L=1, 3, 2 
pO5 J=1, 2 
K=L+J 
KA= 7 - K 
A(K) = AC(KA) + ACK)*W*Z 
W=W+t1.0 


APPRXU = A(2)/A(3) 
APPRXL = A(5)/A(4) 
C = APPRXU - APPRXL 


IF (C.GE.TERMIN) GO TO 4 


x1 
EO 


APPRXL/XO 
=A lhe 2PI“ERP( -035-Y) 


Get accurate value of deviate by Taylor Series. 
(Xl, X2, X3 are derivatives for the Taylor series) 


7 


= XO*=X1*X1 - X1 
= X1**3 + 2. OFXO*XI*XKZ - X2 


= F + LOG(PO) 
= ((X3*D/3. 0+X2)*D/2. O+X1)*D + XO 


D 
X 
X 
X 


PE ( Pee 075)) GO TO100 
= 


100 RETURN 


END 
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APPENDIX D. TABULATED RUN RESULTS FOR RETP 
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Table 1A : 8 Exp in Series, RS = 0.9305 (H1) 
min 4 = 0.0002 f/hr, max A = 0.0016 f/hr, UT = 10 hrs. 


Test Deg. of 
S/N 
Plan Freedom RSLOW LEVEL 
Test 10 until 2NFC 0.9262 0.9630 


ee (160) 0.9256 0.9220 
NFC= 80 2(NFC + 0.9196 0.9970 


NCOMP 
(176) ) 0.9188 0.9910 
2NFC- 0.9296 0.9209 


NCOMP 
(152) 0.9291 0.8430 
2(NFC- 0.9329 0.8400 


NCOMP 
om 0.9325 0.7300 











2 Test 15 until 2NFC 0.9277 0.9550 
15 failed (240) 0.9273 0.9080 
NEC = 120 NEG + 0.9233 0.9899 
(256) 0.9228 0.9750 
2NFC- 0.9298 0.9159 
NCOMP 
(232) 0.9296 0.8329 
2(NFC- 0.9321 0.8440 
NCOMP 
(224) 0.9318 0.7470 
3 Test 15 until 2NFC 0.9268 0.9550 


I1 failed 


(176) 0.9262 0.9159 
2(NFC + 0.9208 0.9960 


NCOMP 
(192) 0.9200 0.9880 
2NFC- 0.9298 0.9159 


NCOMP 
(168) 0.9293 0.8430 
2(NFC- 0.9328 0.8430 


NCOMP 
(160) ) 0.9324 0.7350 


NEC — 55 
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Table 1A : 8 Exp in Series, RS = 0.9305 (Hi) (Cont...) 
min A = 0.0002 f/hr, max A = 0.0016 f/hr, UT = 10 hrs. 


Test Deg. of 
S/N 
Plan Freedom RSLOW LEVEL 
Test 15 until INEC 0.9241 0.9700 


7 failed (112) 0.9229 0.9310 
NFC = 56 2(NFC + 0.9145 0.9980 


NCOMP 
ee 0.9130 0.9940 | - 
2NFC- 0.9289 0.9190 


NCOMP 
= 0920 | 0.8830 
20NEC. | | go “Oe3e2 0.8350 


Test 15 until INFC 
3 dal (48) 
nEc=m | ONREe 


NCOMP 
64) 0.8877 1.0000 
INF "| TOM “a7z68 0.9439 


NCOMP 
- as260 | 0.8600 
2NFC- | 0.1 | 0.9394 0.7530 
NCOMP 
oe assea | 0.6339 





86 


Table 1B : 8 Exp in Series, RS = 0.8025 (Low) 
min 4 = 0.0010 f/hr, max 4 = 0.0045 f/hr, UT = 10 hrs. 


Tes Deftof 
S/N 
Plan Freedom RSLOW LEVEL 
Test 10 until ee fOr | 017907 0.9710 


im tele (160) 0.7884 0.9310 
NFC= 80 ANFC+ | 0.1 | 0.7735 0.9990 


NCOMP) 
(176) 0.7707 0.9920 
2NFC- | 0.1 | 0.7994 0.9299 


NCOMP 
(152) 0.7975 0.8500 
2(NFC- 0.8083 0.8480 


NCOMP) 
(144) 0.8066 0.7410 
INFC | 0.1 | 0.7940 0.9620 
(240) 0.7933 0.9159 

ANFC+ | 0.1 | 0.7826 0.9930 


NCOMP 

(256) ) 0.7816 0.9809 
2NFC- | 0.1 {| 0.7998 0.9240 
NCOMP 

(232) 0.7992 0.8400 
2NFC- | 0.1 | 0.8057 0.8580 
NCOMP 
woo 0.8052 0.7550 
3 Test 15 until 2NFC | 0.1 | 0.7919 0.9660 


Mitayec (176) 0.7896 0.9209 
NFC = 88 ANFC+ | 0.1 | 0.7764 0.9960 


NCOMP 
(192) 0.7735 0.9890 
2NFC- 0.7999 0.9220 


NCOMP 
(168) 0.7978 0.8550 
2NFC- | 0.1 | 0.8079 0.8550 


NCOMP 
(160) ) 0.8061 0.7460 








Test 15 until 
15 failed 






NFC= 120 
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Table 1B : 8 Exp in,Semes, RS =s0ss02>s( Po. Conwe) 
min A = 0.0002 f/hr, inax 4 = 0.0016 f/hr, UT = 10 hrs. 


Test Deg. of 
S/N 
Plan Freedom RSLOW LEVEL 
Test 15 until NEC 0.7851 0.9740 


7 failed (112) 0.7816 0.9380 
2(NFC + 0.7605 0.9980 


NCOMP 
(128) ) 0.7559 0.9960 
2NFC- 0.7978 0.9250 


NCOMP 
(104) 0.7948 0.8600 


2(NFC- 0.8107 0.8410 


NCOMP 
(96) 0.8083 0.7270 
INFC 0.7593 0.9890 
(48) 0.7555 0.9770 
2~NFC+ | 0.1 | 0.7001 1.0000 
NCOMP 
(64) ) 0.6928 1.0000 


NEC] | TOP Giei2 0.9489 


NCOMP 
os 07892) 08720 
ANFC- | 01 | 0.8249 0.7630 


NCOMP 
- oszas | 0.6820 




















NFC= 56 

























Test 15 until 
3 failed 








NFC= 24 
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Table 2A : 8 Wei in Series, RS = 0.9798 (H1) 
min A = 0.001 f/hr, max A = 0.008 f/hr, UT = 10 hrs. 


et Bey. of 
Plan Freedom RSLOW LEVEL 
Test 10 until INFC | 0.1 | 0.9602 0.9859 
10 failed (160) 0.9506 0.9820 
ANFC+ | 0.1 | 0.9565 0.9890 
NCOMP 
176) 0.9461 0.9870 
2NFC- | 0.1 {| 0.9619 0.9820 
NCOMP 
(152) 0.9529 0.9809 
2(NFC- 0.9638 0.9800 
NCOMP 
(144) 0.9553 0.9770 
2 Test 15 until NEC | Ot | 0.9706 0.9770 


IS-dailed (240) 0.9649 0.9719 
NFC= 120 ANFC+ | 0.1 | 0.9687 0.9820 


NCOMP) 
(256) 0.9627 0.9790 
2NFC- | 0.1 | 0.9715 0.9730 

NCOMP 
(232) 0.9659 0.9640 
2ANFC- | 0.1 | 0.9724 0.9650 

NCOMP) 
(224) 0.9671 0.9590 
INFC | 0.1 | 0.9704 0.9719 


(176) 0.9637 0.9650 
A™NFC+ | 0.1 | 0.9679 0.9800 


NCOMP 
-” 0.9606 0.9760 
rie ou i,  O0S7i6 0.9660 

NCOMP 
(168) 0.9652 0.9590 
2(NFC- 0.9729 0.9590 

NCOMP 
(160) 0.9668 0.9510 






NFC= 80 









Test 15 until 
11 failed 






NFC= 88 
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Table 2A : 8 Wei in SeriesaRS — O27 2s Cont.) 
min 2 = 0.001 f/hr, max 4 = 0.008 f/hr, UT = 10 hrs. 
Test Deg. of |___ Measures of accuracy 
S/N 
Plan Freedom RSLOW LEVEL 
Test 15 until 2NFC [01 | 0.9764 | 0.9400 
—_ om 0.9668 
A~NFC+ | 01 | 0.9733 | 0.9579 


NCOMP 
a 0.9624 0.9510 
2NFC- 0.9779 0.9240 
NCOMP 
2(NFC- 0.9795 0.9080 
NCOMP 
0.9713 0.8900 

















NFC= 56 
























(96) 
2NFC 0.9889 0.8000 


(48) 0.9814 0.7770 
2(NFC+ 0.9857 0.8469 


NCOMP 
(64) ) 0.9758 0.8360 
2NFC- 0.9906 0.7510 


NCOMP 
(40) 0.9843 0.7280 
2(NFC- 0.9922 0.6890 


NCOMP 
° 0.9872 0.6530 







Test 15 until 
3 failed 













NEC =24 
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Table 2B : 8 Wei in Series, RS = 0.8323 (Low) 
min 2 = 0.003 f/hr, max 2 = 0.024 f/hr, UT = 10 hrs. 


sy] ies Deg. of 
I Plan Freedom RSLOW LEVEL 


Test 10 until INFC 
10 failed (160) 


NFC= 80 PA es! Oy 
NCOMP) 
(176) 


2NFC- 
NCOMP 
(152) 


2(NEC- 
NCOMP) 
(144) 


y) Test 15 until INFC 
15 failed (240) 


NFC = 120 2NFC+ 
NCOMP) 
(256) 


2NFC- 
NCOMP 
(232) 


2(NEC- 
NCOMP) 
(224) 


5 Test 15 until INFC 
11 failed (176) 


NEC = 88 ZUNE + 
NCOMP) 
(192) 


2NFC- 
NCOMP 
(168) 
2(NFC- 
NCOMP) 
(160) 





91 


a Re iG 0.9800 
0.7383 0.9740 
Per 0.7425 0.9880 


anim | asa 
0.7714 0.9730 


0.7491 0.9680 
0.7813 0.9680 


0.7601 0.9620 
0.8002 0.9640 
0.7766 0.9529 
wi | Oe 0.9730 


a6 | 0.9608 
| O1 | 0.9059 0.9560 


0.7829 0.9430 
"po... o8iic 0.9450 


0.7893 0.9320 
| 0.1 | 0.8056 0.9529 
0.7803 0.9430 
| 0.1 | 0.7909 0.9710 


0.7636 0.9640 
Pom. Weis 0.9430 


0.7888 0.9270 
0.8206 0.9270 


0.7994 0.9060 


Table 2B : 8 Wei in Semess KS = see5z5 } (Low) (Cont...) 
min A = 0.003 f/hr, max 2 = 0.024 f/hr, UT = 10 hrs. 


ss] rm | ess | + aston fever 
: men Eng 
Test 15 until 2NFC 

“——— of 
NFC= 56 2(NFC + 


NCOMP 
| (128) 0.7832 0.9320 
2NFC- 0.8399 0.8800 


NCOMP 
(104) 0.8183 0.8540 
2(NFC- 0.8504 0.8480 


NCOMP 

(96) 0.8304 0.8070 
Die all 0.8801 0.7449 
SL 0.8569 0.7060 

2(NFC+ 0.8475 0.8660 


NCOMP 
(64) 0.8171 0.8410 
2NFC- 0.8970 0.6540 


NCOMP 
40) 0.8778 0.5969 
2(NFC- 0.9145 0.5270 


NCOMP 
at 0.8994 0.4470 


















Test 15 until 
3 failed 






NFC= 24 
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Table 3A : 4 Exp and 4 Wei (mixed) in Series, RS = 0.9801 (Hi) 
min 4 = 0.0002 f/hr, max 2 = 0.0008 f/hr, UT = 10 hrs. 


Test Deg. of 
S/N 
Plan Freedom RSLOW LEVEL 
Test 10 until INFC | 0.1 | 0.9796 0.9190 


10 failed (160) 0.979] 0.8720 
NFC = 80 2NFC+ | 0.1 | 0.9777 0.9740 


NCOMP 
(176) 0.9771 0.9489 
ies | Ol |) 0.9805 0.8720 


NCOMP 
(152) 0.9801 0.7980 
2(NFC- 0.9815 "0.7950 


aun 
2 | Test 15 until NFO OT Os802 | 018950 
oe ey 

NFC= 120 2(NFC + 


NCOMP 
| (256) 0.9786 0.9209 
recy | Ona Oosos 0.8370 


NCOMP 
” (239) 0.9805 0.7530 
ANFC- | 0.1 | 0.9814 0.7670 


NCOMP) 

3 Test 15 until INFC |_O1_ | 0.9802 | 0.8920 
— 
NFC= 88 Cl ht itr air i a 


NCOMP) 
(192) 0.9779 0.9320 
2NFC- 0.9810 0.8410 


NCOMP 
(168) 0.9806 0.7550 
2NFC- | 01 | 0.9819 0.7580 


NCOMP 
(160) 0.9815 0.6590 
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Table 3A : 4 Exp and 4 Wei in Series, RS = 0.9801 (Hi) (Cont...) 
min 4 = 0.0002 f/hr, max 4 = 0.0008 f/hr, UT = 10 hrs. 


Test Deg. of 
Plan Freedom RSLOW LEVEL 
Test 15 until 2NFC 0.9803 0.8950 


tee (112) 0.9797 0.8260 
NFC = 56 2(NFC + 0.9777 0.9690 


NCOMP 
(128) 0.9770 0.9439 
2NFC- 0.9815 0.8060 

NCOMP 
(104) 0.9811 0.7300 
2(NFC- 0.9828 0.7090 

NCOMP 
(96) 0.9825 0.6000 
2NFC 0.9815 0.8730 


(48) 0.9800 0.8020 
2(NFC + 0.9757 0.9740 


NCOMP 
(64) 0.9739 0.9560 
2NFC- 0.9839 0.7250 


NCOMP 
(40) 0.9831 0.6260 
2ANFC- | 0.1 | 0.9868 0.5150 


NCOMP 
(32) ) 0.9862 0.4140 










1 

















Test 15 until 
3 failed 














NFC=24 
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Table 3B : 4 Exp and 4 Wei (mixed) in Series, RS = 0.8089 (Low) 
min A = 0.002 f/hr, max A = 0.008 f/hr, UT = 10 hrs. 


Test Deg. of 
S/N 
Plan Freedom RSLOW LEVEL 
Test 10 until INFC | 0.1 | 0.7957 0.9570 


dln (160) 0.7914 0.9240 
NEC = 80 ANFC+ | 0.1 | 0.7789 0.9920 


NCOMP 
(176) ) 0.7738 0.9770 
2NFC- 0.8042 0.9240 


NCOMP 
(152) 0.8003 0.8730 
2(NFC- 0.8129 0.8710 


Nn 
2 | Test 15 until Now | POUL 08033 is 0.9340 | 
mage ey 

NFC= 120 aNFC+ 


NCOMP 
(736) ) 0.7893 0.9500 
2NFC- 0.8089 0.9010 

NCOMP 
(30) 0.8064 0.8260 
2(NFC- 0.8145 0.8400 

NCOMP 
(224) 0.8122 0.7750 
3 Test 15 until 2NFC — 0.8045 0.9290 


I] failed (176) 0.7993 0.8810 
NFC= 88 2(NFC + 0.7897 0.9850 


NCOMP 
a” 0.7839 0.9570 
2NFC- | 0.1 | 0.8119 0.8820 


NCOMP 
(168) 0.8072 0.8189 
ANFC- | 01 | 0.8195 0.8200 


NCOMP 
(160) ) 0.815] 0.7300 





ws 


Table 3B : 4 Exp and 4 Wei in Series, RS = 0.8089 (Low) (Cont...) 
min A = 0.002 f/hr, max A = 0.008 f/hr, UT = 10 hrs. 


Test Deg. of 
S/N 
Plan Freedom RSLOW LEVEL 
Test 15 until INFC | 0.1 | 0.8056 0.9140 


J fells (112) 0.8017 0.8600 
NFC = 56 A™NFC+ | 0.1 | 0.7830 0.9760 


NCOMP 
C128) 0.7780 0.9610 


2NFC- | 0.1 | 0.8172 0.8460 
NCOMP 
(104) 0.8139 0.7630 


2NFE? O01) SiGraea0 0.7410 


NCOMP 
(96) 0.8262 0.6490 
Fe (er 0.8670 
(48) 0.8085 0.8060 
A~NFC+ | 0.1 | 0.7703 0.9780 
NCOMP. 
(64) ) 0.7571 0.9550 


2NFC- | 0.1 | 0.8424 0.7219 
NCOMP 
(40) 0.8357 0.6050 
ANFC- | 0.1 | 0.8685 0.4859 
NCOMP 
— 0.8641 0.3810 































Test 15 until 
3 failed 





NFC= 24 
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Table 4A : 8 Exponential in Parallel, RS 
min 4 = 0.1000 f/hr, max 4 
Degrees of Freedom = 2NFC 


S/N Test Deg. of 
Plan Freedom RSLOW LEVEL 


Test 10 until 0.9306 0.9290 
10 failed ies 
0.2 0.9320 0.8180 
NFC = 80 


0.9345 (Hi) 
0.1525 f/hr, UT = 10 hrs. 





















Test 15 until 0.9325 0.9190 
15 failed 
240 
0.2 0.9324 0.8210 
NFC = 120 








3 Test 15 until 0.9309 0.9240 
11 failed ae 
0.2 0.9312 0.8329 
NFC= 88 
4 Test 15 until 0.9307 0.9209 
7 failed He 
0.2 0.9308 0.8310 
NFC = 56 
5 Test 15 until 0.9280 0.9220 
3 failed 4g 
0.2 0.9286 0.8370 
NEC = 24 
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Table 4B : 8 Exponential in Parallel, RS = 0.8262 (Low) 
min A = 0.1000 f/hr, max A = 0.2575 f/hr, UT = 10 hrs. 
Degrees of Freedom = 2NFC 


SN] Tes Deg, of 
Plan Freedom RSLOW LEE 


Test 10 until 0.8309 0.8870 


10 failed 
0.2 0.8288 0.7880 
NFC= 80 


Test 15 until 0.8290 0.8870 


15 failed 
0.2 0.8283 0.7860 
NFC= 120 


Test 15 until 0.8266 0.8990 


11 failed 
0.2 0.8264 0.7990 
NFC= 88 


Test 15 until 0.8300 0.8860 


7 failed 
0.2 0.8284 0.7890 
NFC= 56 


Test 15 until 0.8354 0.8810 


3 failed 
NEC=24 0.8364 0.7700 
{ ——4 





98 





Table 5A : 8 Wei in Parallel , RS = 0.9265 (Hi) 
min A = 0.100 f/hr, max A = 0.128 f/hr, UT = 10 hrs. 
Degrees of Freedom = 2NFC + 0.5(NC/NF) 


S/N Test Deg. of 
— Freedom RSLOW LEVEL 
Test 10 until 0.9154 0.9540 
10 failed 16] 


0:97 0.8880 


Test 15 until 0.9179 0.9510 


15 failed 7A] 
0.2 0.9170 0.8770 
NFC= 120 


Test 15 until 0.9248 0.9180 






















3 
11 failed 177 
0.2 0.9250 0.8130 
NFC= 88 
Test 15 until 0.9323 0.8630 
7 failed 113 
0.2 0.9318 0.7589 
NFC = 56 






0.9152 0.9270 


3 
0.8980 0.8850 


Test 15 until 
3 failed 





oN 


Table 5B : 8 Wei in Parallel , RS = 0.8351 (Low) 
min 2 = 0.100 f/hr, max 4 = 0.163 f/hr, UT = 10 hrs. 
Degrees of Freedom = 2NFC + 0.5(NC/NF) 


S/N Test Deg. of 
Plan Freedom RSLOW Lauer 













Test 10 until 0.8296 0.9180 
10 failed 161 
0.2 0.8331 0.8140 
NFC = 80 


Test 15 until 0.8322 0.9100 | 









15 failed A} 
0.2 0.8314 0.8279 
NFC= 120 
Test 15 until 0.8494 0.8380 
11 failed 177 
0.2 0.8489 0.7130 
NFC= 88 


Test 15 until or | ers 0.7560 


7 failed ve 
0.2 0.8705 0.6450 
NFC= 56 
Test 15 until 0.8678 0.8360 
3 failed 5] 


0.8469 0.7770 








NFC= 24 
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Table 6A : 4 EXP and 4 Wei (mixed) in Parallel, RS = 0.9408 (H1) 
min A = 0.100 f/hr, max 2 = 0.130 f/hr, UT = 10 hrs. 
Degrees of Freedom = 2NFC + 0.5(NC/NF) 


S/N Test Deg. of 
Plan Freedom RSLOW LEW 


Test 10 until 0.9342 0.9420 


10 failed 
0.2 0.9366 0.8500 
NFC= 80 


Test 15 until 0.9362 0.9330 




















15 failed Bs 
0.2 0.9363 0.8510 
NFC= 120 
Test 15 until | 0.1 | 0.9389 0.9240 
11 failed 177 
0.2 0.9378 0.8310 
NFC = 88 


Test 15 until ce ic 0.8920 


7 failed 113 
02 0.9427 0.7819 
NFC= 56 


Test 15 until Con ores 0.9080 


5 
3 failed 
0.2 0.9325 0.8410 
NFC= 24 
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Table 6B : 4 EXP and 4 Wei (mixed) in Parallel, RS = 0.8170 (Low) 
min A = 0.100 f/hr, max 4. = 0.220 f/hr, UT = 10 hrs. 
Degrees of Freedom = 2NFC + 0.5(NC/NF) 


RSLOW LEVEL 


Test 10 until 
10 failed 


IN cm 
Test 15 until 20.1 ol... Ofa62 0.8630 


15 failed 
0.2 0.8236 0.7560 
NFC = 120 


Test 15 until Te os: 0.8290 


11 failed 
0.2 0.8334 0.7150 
NFC= 88 


Test 15 until on) ees: 0.7660 


7 failed 
0.2 0.8469 0.6470 
NFC= 56 


Test 15 until | 0.1 | 0.8650 0.7819 


3 failed 
0.8563 0.6709 
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Table 7A : 8 Exp in Series- Parallel, RS = 0.9249 (Hi) 
min 2 = 0.0003 f/hr, max 4 = 0.0024 f/hr, UT = 10 hrs. 
Degrees of Freedom = 2NFC 


S/N Test Deg. of 
Plan Freedom RSLOW BeweEL 


Test 10 until 
10 failed 


Test 15 until 
15 failed 
NFC= 120 
Test 15 until 
11 failed 
NEFC = 88 
Test 15 until 
7 failed 
NFC= 56 


Test 15 until 
3 failed 


160 


240 


176 


112 


for | "oral 0.9159 


0.9226 0.8620 


0.9253 0.8960 


0.9236 0.8360 


0.9252 0.8950 


osms | 0.8390 


0.9222 0.9260 


0.2 0.9204 0.8789 


ee ey 0.9620 


09123 | 0.9260 
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Table 7B : 8 Exp in Series- Parallel, RS = 0.8228 (Low) 
min 4 =0.00075 f/hr, max 2 =0.0060 f/hr, UT =10 hrs. 
Degrees of Freedom = 2NFC 


SN] Test a 
Plan Freedom RSLOW LEVEL 
Test 10 until 0.8210 0.9159 


10 failed 
0.2 0.8175 0.8620 
NFC= 80 


Test 15 until 0.8237 0.8960 


15 failed 40 
0.2 0.8199 0.8360 
NFC = 120 
= Test 15 until | 0.1 | 0.8234 0.8950 


















11 failed 176 
0.2 0.8176 0.8390 
NFC= 88 


Test 15 until | 0.1 | 0.8168 0.9260 


7 failed 112 | 
0.2 0.8127 0.8789 
NEC =56 


Test 15 until | 0.1 | 0.8015 0.9620 


3 failed 
0.2 0.7949 0.9260 
NFC=24 | 
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Table 8A : 8 Wel in Series- Parallel, RS = 0.9328 (H1) 
min A = 0.002 f/hr, max A = 0.016 f/hr, UT = 10 hrs. 
Degrees of Freedom = 2NFC + 0.5(NC/NF) 


Test Deg. of 
Plan Freedom RSLOW LEVEL 
Test 10 until boi || Side 0.9579 


10 failed 
_ os 0.8939 0.9510 


Test 15 until | 0.1 | 0.9247 0.9330 


15 failed 
0.2 0.9108 0.9230 
NFC= 120 


Test 15 until | 0.1 | 0.9266 0.9220 


11 failed 
0.2 0.9149 
NFC= 88 


Test 15 until | O04 | 0.9414 0.8609 


7 failed 
0.2 0.9250 0.8419 
NEFC= 56 


Test 15 until | 0.1 | 0.9699 0.7079 


3 failed 
0.2 0.9536 0.6740 
NFC = 24 
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Table 8B : 8 Wel in Series- Parallel, RS = 0.8321 (Low) 
min A = 0.00325 f/hr, max A = 0.026 f/hr, UT = 10 hrs. 
Degrees of Freedom = 2NFC + 0.5(NC/NF) 


Test ener 
Plan Freedom RSLOW LEVEL 
Test 10 until | 0.1 | 0.7963 0.9560 
10 failed 


0.2 0.7703 0.9470 
NFC = 80 
Test 15 until 0.8204 0.9310 


15 failed 7A) 
0.2 0.7986 0.9200 
NFC= 120 






















3 Test 15 until 0.1 | TORREG 0.9190 
11 failed Ga 
0.2 0.8046 0.9000 
NFC= 88 
Test 15 until 0.8485 0.8590 
7 failed 113 
0.2 0.8237 0.8329 
NFC= 56 









| 0.1 | 0.8999 0.7500 


Sl 
0.8599 0.6950 


Test 15 until 
3 failed 
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Table 9A : 4 EXP and 4 Wei in Series- Parallel, RS = 0.9276 (Hi) 
min A = 0.005 f/hr, max 4 = 0.020 f/hr, UT = 10 hrs. 
Degrees of Freedom = 2NFC + 0.5(NC/NF) 


S/N Test Deg. of 
Plan Freedom RSLOW LEVEL 
Test 10 until / 0.1 | 0.9096 0.9400 


10 failed 
0.2 0.8904 0.9349 
NFC= 80 


Test 15 until Pou” jl diiszes 0.9170 


















15 failed ae 
0.2 0.9079 0.9040 
NFC= 120 
Test 15 until | 0.1 | 0.9266 0.9050 
11 failed 177 
0.2 0.9109 0.8850 
NFC = 88 







Test 15 until Poi | Oe 0.8550 
7 failed 113 
0.2 0.9203 0.8390 
NFC= 56 


Test 15 until | 0.1 | 0.9637 0.7430 


5 
3 failed 
0.2 0.9451 0.7200 
NFC= 24 
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Table 9B : 4 EXP and 4 Wei in Series- Parallel, RS = 0.8248 (Low) 
min 2 = 0.008 f/hr, max 4 = 0.032 f/hr, UT = 10 hrs. 
Degrees of Freedom = 2NFC + 0.5(NC/NF) 


S!N Test eower 
Plan Freedom RSLOW LEVEL 


Test 10 until | 0.1 | 0.7909 0.9420 


10 failed 161 
0.2 0.7627 0.9360 
NFC= 80 


Test 15 until | 0.1 | 0.8164 0.9150 


15 failed 
241 
0.2 0.7915 0.9020 
NFC= 120 


Test 15 until Or | Waa 0.9010 


11 failed 177 
0.2 0.7964 0.8820 
NFC= 88 























Test 15 until "On. | eRe 0.8660 
7 failed 113 
0.2 0.8119 0.8419 
NFC= 56 









| 0.1 | 0.8683 0.8170 


Test 15 until 
3 failed 


NFC= 24 
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Table 10A :10 EXP and 5 Wei in Series- Parallel, RS = 0.9472 (H1) 
Exp: min 2 = 0.025 f/hr, max A = 0.075 f/hr. 
Wei: min? = 0.055 f/hr, max 4 = 0.095 f/hr. 
UT = 10 hrs. 
Degrees of Freedom = 2NFC + 0.5(NC/NF) 


S/N Test Deg. of 
viet Freedom RSLOW LEVEL 
Test 10 until 0.9383 0.9710 


10 failed 30] 
0,2 0.9369 0.9550 
NFC= 150 


Test 20 until | 0.1 | 0.9449 0.9370 


20 failed 
0.2 0.9432 0.8950 


Test 20 until | 0.1 | 0.9444 0.9349 










NFC= 300 


15 failed 45] 
: 0.2 0.9430 0.8979 
NFC= 225 


Test 20 until | 0.4 | 0.9435 0.9290 


10 failed 30] 
0.2 0.9415 0.9000 
NFC= 150 


| 0.1 | 0.9344 0.9809 


152 
0.9310 0.9520 





Test 20 until 
5 failed 






NFC=75 
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Table 10B :10 EXP and 5 Wei in Series- Parallel, RS = 0.8324 (Low) 
Exp : min * = 01025%/hr; maext =—"OM7S Thr. 
Wei: min 4 = 0.055 f/hr, max 2 = 0.125 f/hr. 
UT = 10 hrs. 
Degrees of Freedom = 2NFC + 0.5(NC/NF) 
ee 
esi ere 
Test 10 until 0.1 | 0.8103 | 09665 


10 failed 30] 
0.2 0.8081 0.9430 
NFC= 150 


Test 20 until 0.8247 0.9380 


20 failed 
0.2 0.8219 0.8990 
NEC = 300 


Test 20 until | 0.1 | 0.8236 0.9470 


15 failed 451 ~ 
0.2 0.8204 0.8990 
NFC= 225 


Test 20 until 0.8232 0.9500 




















10 failed 30] 
0.2 0.8166 0.8950 
NFC= 150 
Test 20 until | 0.1 | 0.8008 0.9760 


5 failed 
0.2 0.7917 0.9489 
NFC=75 
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APPENDIX E. TABLE OF CHI-SQUARE DISTRIBUTION 


Table 21 in this appendix provides the eightieth and ninetieth 
percentile points for the chi-square distribution with degrees of freedom 
ranging from | to 402 in increments of one. The percentile points appear 
under the column headed CHI in Table 21. 

The computational algorithm used to construct these percentile points 
is defined in the Chi-Square Quantile Function routine together with the 
Standard Normal Variate Computation routine on pages 79 - 83 in this 
thesis. 

Values of these chi-square percentile points for degrees of freedom as 
large as 601 were needed in this thesis. The extensive set of chi-square 
percentile values in Table 21 is provided as a convenience for a reader 


who may want to apply this methodology to a particular system. 
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Table 21 : Chi-square distribution 
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Table 21 : Chi-square distribution (Continued...) 


Cc egeeeamweeee@e@eGgeeeeeetewaeewe eee es eweww et eee2 ete et eB ew Bw BZwBVZV—Ss SBT BZ FBZ FZVWSs SB HBs BZ Bsw BZ BZ FZT ZT SBVTSsT STS SFT TBTZVW_sZT_s @ = 
Left tail area = 0.9 Left tail area = 0.8 
[a =Eueueueueeue @Seeweeees Ss Besse eweswesw sews es Bestest es FB Bes ess FB FBI SFB FBT BWI es SFT FT TIE FTBwwBVw—_sswZwsBZw_estQ7° gg sBeee@t FB FFB SB FI BBs & 
DF CHI DF CHI DF CHI DF CHI DF CHI DF CHI 
Cc [ee weeeweuztew @&eeew ee eee eee ete eet estes SFB FBT SFB FB FBT FBI FBI sBI FT FT FI TFT TFT Bee qg@sweerae@qe@eegw@er@#s8w=—eBgeqg@esrstsws#uesea = =s 
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Table 21 : Chi-square distribution (Continued...) 
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